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ABSTRACT 

A  laser- source   integrating   sphere   reflectometer   was   designed   and  built 
to   measure    the    reflectance    of   specimens    at  high  temperatures.       The    sphere 
was    calibrated   for   linearity  of  response   at  0.  632y.  by  means    of  a   shallow 
cylindrical   cavity  with  a  variable   depth- to- radius   ratio,    having   a  lining   of 
known  reflectance.       Preliminary  tests    showed  that  the    flux   emitted  by   a  hot 
specimen  at  temperatures   up  to   2500    "K  will  not   invalidate  the   reflectance 
measurement.      An     ellipsoidal  mirror   reflectometer   was   calibrated   for   all 
known  errors    in  the   1  to   7.  5)i  range.       A  preliminary   analysis    indicates   that 
the   errors   in  measurement   of  absolute   reflectance   with  this    equipment    should 
not  exceed  2   percent.       A  review  of  the   literature   on  relation  of  thermal 
radiation  properties   to   other   properties    of  materials   is   presented,    together 
with  a   summary   of  the   work  done   in  an  effort  to   compute   the    spectral 
emittance   of  rhodium. 

Key  words:      Emissivity,    emittance,    high  temperature   reflectance,    infrared 
reflectance,    radiation  properties,    reflectance,     spectral   emit- 
tance,   spectral  reflectance,    thermal  radiation. 
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Symbols  (Cont'd) 

e(cp,A.)  spectral  emissivity  in  direction,  cp  at  wavelength  X 

6  phase  component  of  complex  quantity,  V 

X  wavelength,  or  wavelength  of  incident  wave 

*o  cAlo  (^)2 

\f  2ncT  ,  viscous  damping  parameter  for  free  electrons 

Xf  viscous  damping  parameter  for  free  electrons,  type  g 

A^h  viscous  damping  parameter  for  bound  electrons,  type  h 

X  ,  resonant  wavelength  of  bound  electrons,  type  h,  with  no  damping  forces 
present 

M-  microns 

M>*  magnetic  permeability 

"  3.14159  ... 

P  reflectance 

Pc  reflectance  of  a  cavity 

Pm  resistivity  of  a  metal 

p  reflectance  of  specimen 

pg  reflectance  of  sphere  wall 

p(0,\)  specular  reflectance  for  normal  incidence,  at  wavelengt-r.  } 

p(cp,A)  specular  spectral  reflectance  for  radiation  incident  at  angle  cp 

pg-t  reflectance  of  standard 

Pw  reflectance  of  cavity  wall 

PWjtpj^i!  reflected  component  polarized  normal  to  the  plane  of  incidence 

p(i|r,cp,A;p  reflected  component  polarized  parallel  to  the  plane  of  incidence 

T  relaxation  time 

9  angle  of  refraction,  measured  from  the  normal 

*  angle  of  incidence,  measured  from  the  normal 

w  frequency 
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1.  Objectives 

The  over-all  objective  of  this  project  was  to  develop  equipment  and  procedures  for 
obtaining  reliable  measurements  of  the  thermal  radiation  properties  of  materials,  and  to 
select,  prepare,  and  calibrate  physical  standards  for  use  by  other  laboratories  in  checking 
their  equipment  and  procedures.   The  specific  temperature  range  of  interest  was  800  to  2500  °K 
(about  1000  to  4000  °F).   Procedures  applicable  to  both  metallic  and  non-metallic  (ceramic) 
materials  heated  in  air,  controlled  atmosphere  and  vacuum  were  desired. 

The  specific  objectives  for  the  report  period  included  the  following,  (l)  A  study  of  the 
feasibility  of  evaluating  normal  spectral  emittance  at  high  temperature  by  measurement  of 
reflectance  by  means  of  a  laser-source  integrating  sphere  reflectometer.  Reflectance  will 
be  measured  at  a  few  specific  wavelengths,  and  values  at  intermediate  wavelengths  will  be 
interpolated  from  the  continuous  spectral  reflectance  curve  of  the  material  measured  at  room 
temperature  with  a  different  instrument.   (2)  Procurement,  preparation,  and  calibration  of 
standards  of  normal  spectral  emittance  in  the  temperature  range  of  1400  to  1900  °K  (approxi- 
mately 2000  to  3000  °F).   (3)  Extension  of  the  wavelength  range  of  the  normal  spectral  emit- 
tance equipment,  in  the  800  to  1400  °K  (approximately  1000  to  2000  °F)  temperature  range  to 
wavelengths  of  15  to  35  microns. 

Development  of  equipment  and  procedures  for  measurement  of  normal  spectral  emittance  by 
the  rotating  cylinder  method  and  of  total  normal  emittance  by  the  shallow  cavity  method  has 
proceeded  on  another  project  sponsored  by  the  George  C.  Marshall  Space  Flight  Center  of  NASA. 
Reports  covering  work  on  both  projects  are  sent  to  both  sponsors. 

2.  Laser-Source  Reflectometer 

2.1  Background 

Spectral  emittance  measurements  obtained  by  direct  comparison  of  the  spectral  radiant 
flux  density  of  a  hot  specimen  to  that  of  a  laboratory  blackbody  furnace  at  the  same  tempera- 
ture are  subject  to  large  errors,  not  only  at  wavelengths  that  are  appreciably  shorter  than 
that  at  the  peak  of  the  Planck  distribution  function  at  the  test  "temperature,  (see  Figure  1, 
Reference  56),  but  also  at  temperatures  above  those  at  which  standard  thermocouples  can  be 
used. 

The  errors  at  the  short  wavelengths  arise  from  two  sources:   (l)  The  energy  available 
for  measurement  decreases  rapidly  with  decreasing  wavelength  at  wavelengths  below  the  peak; 
hence,  there  is  not  enough  energy  to  permit  measurements  of  the  highest  accuracy  at  the 
shortest  wavelengths.   (2)  The  temperature  dependence  of  spectral  flux  density  increases 
rapidly  with  decreasing  wavelength  at  wavelengths  below  the  peak;  therefore,  a  small  temper- 
ature difference  between  specimen  and  blackbody  reference  will  cause  a  much  larger  error  in 
emittance  at  wavelengths  shorter  than  the  peak  than  at  wavelengths  longer  than  the  peak* 

At  temperatures  above  those  at  which  standard  thermocouples  can  be  used,  the  temperature 
scale  is  based  upon  the  optical  pyrometer,  and  either  an  optical  pyrometer  or  a  radiation 
pyrometer  is  commonly  used  for  temperature  measurements  in  this  range.  These  instruments 
are  accurate  for  measuring  the  temperature  of  a  laboratory  blackbody  furnace,  and  they  also 
measure  the  brightness  temperature  or  the  radiance  temperature  of  a  specimen  with  acceptable 
accuracy.  Conversion  of  this  brightness  or  radiance  temperature  to  true  temperature,  however, 
(requires  knowledge  of  the  spectral  emittance  of  the  specimen  over  the  appropriate  wavelength 
range.  Thus  the  emittance  must  be  known  in  order  to  measure  the  temperature,  which  in  turn 
must  be  known  to  evaluate  the  emittance. 

Reflected  flux,  as  contrasted  to  emitted  flux,  is  nearly  independent  of  the  temperature  of 
the  specimen,  at  least  over  temperature  intervals  that  do  not  exceed  about  50 °K.  For  opaque 
specimens,  the  emittance  is  equal  to  one  minus  the  reflectance,  by  Kirchhoff 's  law.  The  meas- 
urement of  reflectance  involves  only  the  comparison  of  incident  and  reflected  flux.  Thus,  such 
measurements  would  appear  to  offer  the  possibility  of  considerably  increasing  the  precision 
with  which  emittance  can  be  determined,  particularly  in  the  short-wavelength  and  high-tempera- 


The  nomenclature  used  in  this  report  is  that  used  in  references  51  through  56,  and  is  defined 
in  detail  in  references  51  and  54. 
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ture  regions  where  direct  measurements  are  subject  to  large  errors.  However,  for  reflectance 
to  be  equal  to  one  minus  the  normal  spectral  emittance,  the  reflectance  must  be  measured  under 
conditions  of  either  normal  illumination  and  hemispherical  viewing,  or  diffuse  illumination 
and  normal  viewing. 

Three  general  methods  have  been  used  to  measure  reflectance  under  these  conditions; 
by  using  (l)  integrating  sphere,  (2)  integrating  hemisphere,  and  (3)  hemispherical  source 
instruments.  Integrating  sphere  instruments  are  generally  preferred,  and  are  the  type  most 
widely  used  at  wavelengths  for  which  suitable  sphere  coatings  are  available.  The  principal 
problem  in  measuring  reflectance  of  a  specimen  at  high  temperatures  is  to  distinguish  between 
the  flux  reflected  by  the  specimen  and  that  emitted  by  the  hot  specimen.  This  can  be  done  by 
chopping  the  beam  of  incident  flux,  so  that  an  a-c  signal  is  produced  by  the  detector,  which 
is  synchronously  amplified  and  recorded.  The  d-c  signal  from  the  flux  emitted  by  the  hot 
specimen  is  not  amplified,  and  appears  only  as  noise.  The  signal-to-noise  ratio  can  be  in- 
creased by  using  a  source  of  high  flux  density  within  the  wavelength  interval  of  interest.  A 
helium-neon  continuous  wave  laser  was  selected  for  this  purpose,  because  its  flux  density 
within  the  narror  wavelength  band  at  which  it  lases  is  several  orders  of  magnitude  larger 
than  the  flux  emitted  by  any  solid  specimen  at  a  temperature  below  its  melting  point,  within 
the  same  wavelength  interval. 

2.2  Theory  of  Integrating  Spheres 

The  ideal  integrating  sphere  for  use  in  reflectance  measurements  consists  of  a  hollow 
sphere  with  three  small  openings  in  the  wall.  The  inside  wall  has  a  completely  diffusing 
surface  of  high  reflectance.  The  specimen  is  placed  over  one  of  the  openings,  called  the  • 
specimen  port,  in  a  position  such  that  it  forms  a  portion  of  the  sphere  wall.  It  is  illu- 
minated by  a  beam  of  light  projected  onto  its  surface  through  a  second  opening,  called  the 
entrance  port.  A  detector  views  a  portion  of  the  sphere  wall  through  the  third  opening, 
called  the  detector  port.  The  flux  initially  reflected  by  the  specimen  is  multiply  reflected 
by  the  sphere  wall  and  uniformly  distributed  over  its  surface.   If  the  detector  and  the  area 
of  the  sphere  wall  that  it  views  are  shaded  from  the  specimen,  the  response  of  the  detector 
will  be  proportional  to  the  reflectance  of  the  specimen,  regardless  of  the  geometric  dis- 
tribution of  the  reflected  flux  (assuming  that  a  negligible  fraction  is  lost  out  the  open- 
ings). 

When  a  standard  of  known  reflectance  is  substituted  for  the  specimen,  the  response  of 
the  detector  is  likewise  proportional  to  the  reflectance  of  the  standard,  and  the  reflectance : 
of  the  specimen  can  be  computed  from  the  known  reflectance  of  the  standard  and  the  measured 
ratio.  This  is  usually  referred  to  as  operation  in  the  substitution  mode. 

Unfortunately,  all  integrating  spheres  depart  from  the  ideal  conditions  to  some  extent, 
and  corrections  must  be  made  for  the  multiply-reflected  flux  lost  out  of  the  openings,  and 
that  absorbed  by  the  specimen.  Such  corrections  can  be  kept  small  by  keeping  the  area  of 
the  openings  small  compared  to  the  total  inside  area  of  the  sphere.  Also,  no  real  sphere 
coating  is  completely  diffusing,  and  the  reflectance  of  the  sphere  wall  may  vary  over  its 
surface.  However,  in  spite  of  rather  large  departures  from  ideal  conditions,  carefully 
constructed  integrating  sphere  reflectometers  do  give  useful  reflectance  data. 

The  theory  of  the  integrating  sphere  reflectometer  for  a  perfectly  diffuse  sphere  wall 
and  diffusely  reflecting  specimens  has  been  presented  by  Jacquez  and  Kuppenheim  [1]=/.  For 
measurements  by  the  substitution  method  (in  which  first  the  standard  and  then  the  specimen 
is  placed  over  the  specimen  port  and  measured)  the  reflectance  of  the  specimen,  P3,is  ex- 
pressed by 
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Figures  in  brackets  refer  to  papers  cited  in  Section  8  References. 


where 

Bs  =  brightness  at  detector  port  with  specimen  at  specimen  port. 

Bg^  —  brightness  at  detector  port  with  standard  at  specimen  port. 

p  =  reflectance  of  sphere  wall. 

Pe  =  spnerical  area  of  entrance  port. 

P(j  =  spherical  area  of  detector  port. 

p  =  effective  spherical  area  of  specimen  port. 

2 
S  =  total  area  of  sphere  wall  IfiR 

S'  =  S  -  pe  -  Pj  -  Ps  total  reflecting  area  of  sphere  wall. 

C  =  calibration  factor  for  sphere,   [~,     st  "  s   s  "1 

L     "   1   -   pS'/s  -   pp  /S     J    ' 
s    s 

I  The  equation  for  C  reduces  to 

i  -  ps'/s  -  Pstps/S 

C  =  1  -  p  s'/s  -  p  p  /s  ^ 

s  rs 

The  ratio  B  /B  +  is  the  ratio  of  the  reflectometer  readings  for  the  specimen  and 

|  standard,  respectively,  and  C  is  the  correction  factor  for  flux  losses  through  the  sphere 

ports. 

The  theory  given  above  is  for  the  perfect  sphere,  where  the  specimen  and  standard  have 
a  spherical  surface  of  the  same  radius  as  the  sphere,  and  are  diffuse.   It  can  be  modified 
for  the  case  of  flat  specimens.  No  theory  is  available  for  specularly  reflecting  specimens. 

2.3  Description  of  Instrument 

The  integrating  sphere  reflectometer  that  was  designed  for  the  present  work  as  described 
in  reference  55  is  shown  schematically  in  Figure  1.  It  consists  of  a  source,  the  integrating 
sphere  proper,  and  a  detector,  plus  auxiliary  equipment. 

A  helium-neon  continuous-wave  gas  laser  is  used  as  the  source.   It  produces  a  slightly 
diverging  beam  of  plane-polarized  coherent  radiation  about  2mm  (0.08  in.)  in  diameter  at  the 
source.  The  laser  will  generate  narrow-band  radiation  at  three  different  wavelengths, 
0.6328,  1.15*  and  3-39  microns.  Because  the  dielectric-coated  mirrors  reflect  the  desired 
wavelength  and  absorb  the  others,  the  lasing  action  is  confined  to  the  desired  single  wave- 
length. Total  power  output  is  about  one  milliwatt.  When  the  power  input  is  properly  ad- 
justed, the  laser  output  in  a  dust-free  atmosphere  is  stable  to  better  than  0.5  percent  for 
periods  of  several  hours.   Input  power,  however,  is  critical.  When  the  power  is  not  pro- 
perly adjusted,  the  laser  output  varies  sinusoidally  with  an  amplitude  variation  of  2  to  5 
percent  and  a  period  that  varies  from  about  10  minutes  to  an  hour  or  more,  depending  on  the 
power  setting. 

The  Brewster-angle  windows  on  the  gas-discharge  tube  collect  a  high  static  charge,  so 
that  dust  paritcles  in  the  atmosphere  are  attracted  to  them.  Such  dust  particles  on  the 
windows  seriously  reduce  the  power  output  of  the  laser.  In  order  to  obtain  the  desired 
stability,  it  was  necessary  to  operate  the  laser  in  a  dust-tight  enclosure. 

The  laser  beam  passes  out  of  the  enclosure  through  a  window,  after  which  it  is  chopped 
at  13  cps.  The  chopped  beam  is  then  directed  downward  through  the  entrance  port  of  the 
sphere  to  strike  the  specimen  at  an  angle  of  about  6°  to  the  normal.  The  beam  is  about  7.5mm- 
in  diameter  (0.3  in.)  at  the  specimen,  but  can  be  reduced  to  a  diameter  of  about  2mm  (0.075 
in.)  by  means  of  a  long  focal  length  lens  inserted  in  the  system  outside  the  sphere. 
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Figure  1.   Schematic  Diagram  of  High-Temperature  Integrating  Sphere 
Reflectometer,  Showing  Auxiliary  Equipment  Required. 
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Figure  3-  Reflectance  of  3M  Velvet  White 
Coating,  Relative  to  MgO. 


Figure  2.  Cross  Sectional  View  of  Inte- 
grating Sphere,  Showing  Location  of 
Ports.   A  shallow  cavity  is  shown  at 
the  specimen  port. 


A  cross-sectional  view  of  the  integrating  sphere  is  shown  in  Figure  2.  The  sphere  coat- 
Lsed  for  the 
vet  white  paint. 


ing  used  for  the  preliminary  measurements  was  a  white  primer  with  a  thin  top  coat  of  3M  Vel- 


The  reflectance  of  the  3M  Velvet  white  sphere  coating  relative  to  freshly  smoked  MgO  is 
shown  in  Figure  3-  Preliminary  goniophotometric  measurements  indicated  that  the  Velvet  coat- 
ing is  a  good  approximation  to  a  perfect  diffuser. 

The  significant  dimensions  of  the  sphere  for  use  with  equation  (2)  are  as  follows: 

pe  =  spherical  area  of  entrance  port,  0.994  in.  (1-1/8  in.  diameter) 

p,  =  spherical  area  of  detector  port  0.442  in.   (3/4  in.  diameter) 

p  =  effective  spherical  area  of  specimen  port  1.227  in.   (1-1/4  in.  diameter) 

S  =  Total  area  of  sphere  wall,  IpR^,   615.752  in  .  Radius  of  sphere  is  7.0  in. 

S!  =  S  -  pe  -  p^  -  ps  total  reflecting  area  of  sphere  wall,  613.089  in. 

p  =  reflectance  of  sphere  wall,  0.825  at  0.6328u. 

Substituting  the  above  values  in  equation  (2),  we  obtain  for  the  sphere,  calibration 
factor  C 

0.1786  -  0.00199pQ+ 

C  =  —  =  1.00 

0.1786  -  0. 00199 p 

where  p  g-^is  the  reflectance  of  the  standard  and  p  s  is  the  reflectance  of  the  specimen.   If 
P_st  and  Pst   are  eq.ual,(-'  =  1-00.  If  pg-t,  =  0.80  and  P  s  =  0,40,  which  is  about  as  large  a  dif- 
ference in  p  s  and  pst  as  is  likely  to  be  encountered,  C  =  0.9955.  The  minimum  value  of  C, 
which  will  never  be  attained,  would  be  0.9889  for  p  ^  =  1.0  and  ps  =  0.  Because  C  deviates 
only  slightly  from  a  value  of  1.00,  it  was  neglected  in  the  preliminary  checking  of  the 
reflect ometer. 

The  detector  used  at  0.6328  micron  was  a  1P28  photomultiplier.   Several  different  detec- 
tor configurations  were  used,  as  are  described  later.  The  signal  from  the  detector  was  am- 
plified and  rectified  by  standard  Perkin-Elmer  electronics  and  recorded  by  a  Leeds  and  North- 
rup  recording  potentiometer.  A  lead  sulphide  detector  and  preamplifier  were  obtained  and  are 
available  for  use  at  1.15  and  3 .39  microns. 

A  low-temperature  heater,  for  heating  specimens  12.7mm  (l/2  in,)  in  diameter  to  tempera- 
tures up  to  about  800  °K  (1000  °F)  has  been  designed  and  built.  A  schematic  drawing  of  the 
heater  is  shown  in  Figure  4.  An  important  design  feature  is  the  provision  of  holders  for 
specimens  of  thicknesses  of  6.5  mm  (1/4  in.)  for  nonmetals  and  1.6  mm  (0.065  in.)  for  metals. 
This  heater  is  for  preliminary  checking  of  the  equipment.  A  heater  for  operation  at  high 
temperatures  (up  to  about  2500  °K  or  4000  °F)  was  not  constructed  during  the  report  period. 

2.4  Calibration  of  the  Integrating  Sphere  Reflect ometer 

2.4a  Linearity  Check 

The  linearity  of  the  response  of  the  detector-amplifier-recorder  system  was  checked  by 
means  of  a  rotating  disc  attenuator  with  blades  that  have  transmittances  of  25.25  +  0.14, 
50.00  +  0.02,  and  75.10  +  0.2  percent.  At  0.6328  micron  it  was  found  necessary  to  insert 
into  the  optical  path  a  filter  having  a  transmittance  of  about  0.005,  in  order  to  avoid  sat- 
uration of  the  detector  and  a  consequent  non-linear  response.  The  filter  was  located  between 
the  integrating  sphere  and  the  detector,  so  that  the  background  radiation  from  a  hot  specimen^ 


•£/  The  primer  was  DuPont  65-3010  white,  sprayed  to  a  thickness  of  about  0.075mm  (0.003 
in.)  and  air  dried. 


List  and  Description  of  Parts 

1  Integrating  sphere  (14"  OD)  &  6"  OD  flange,  flange  outer  flat 
surface  is  approximately  l/2"  below  sphere  at  heater  centerline 

2  Heated  silver  core  with  cavities  for  l/2"  diameter  specimens 
of  .070"  and  .250"  thicknesses,  TC  hole  diameter  .099" 

3,4,5  Three  layers  of  GE  asbestos  paper,  5  ft.  18  gauge  (.050"D) 

nichrome  wire  at  .425  n/ft.,  aluminum  phosphate  cement  coating 

6,  7   Mica  spacers  and  washers,  clearance  at  6  is  .030" 
clearance  at  7  is  3/32" 

8  End  plate,  l/2"  thick,  transite  coat  with  SiOg  -  NagSiOa 

9  Tie  down  bolt,  l/8"  dia.  nominal,  3  required  at  120°,  steel 

10  Sphere  wall  filler  disc,  .065  thick,  aluminum,  coat  with  SiO^ , 
Na-,Si03 

11  Alignment  boss^  aluminum 

12  Cap  screw  10-32,  3  required  at  120°  with  washer 

13  Body  assembly,  brass 

14  Cooling  coils,  1/8"  0D  tubing 

15  Cap  screw  8-32,  3  required  at  120° 

16  Insulation,  Fiberfrax 

17  Lower  cap  and  cover,  brass 


Figure  4.  Specimen  Heater  Designed  to  Heat  12.7  mm  (g  Inch)  Diameter 
Specimens  to  Temperatures  Up  to  About  800  °K  (1000  °F).  Specimens  of 
metals  1.6  mm  (0.065  inch)  thick  can  be  accommodated  with  the  heater 
in  the  position  shown.  By  reversing  the  heater,  dielectric  specimens 
6.35  mm  (£  inch)  thick  can  be  accommodated. 
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Figure  5.   a.  First  Reflection  Cavity.  The  base  of  the  cavity  is 
a  diffuse  reflector,  but  the  walls  are  lined  with  black  velvet, 
so  that  only  the  flux  directly  reflected  by  the  base  is  permitted 
to  escape  from  the  cavity  opening,  b.  Lined  Cavity.  Both  base 
and  walls  of  the  cavity  reflect  diffusely. 


The  moveable  plunger  in  each  cavity  permits  the  depth  to  radius 
ratio  (L/R)  to  be  varied  from  0.0  to  2.0.  3M  Velvet  yellow  paper 
was  used  as  the  diffusely  reflecting  surfaces. 
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Assembly  for  Integrating  Sphere  Reflec- 
tometer. 


Figure  6.  Areas  of  Sphere  Directly 

Illuminated  by  Flux  Reflected  from  First 
Reflection  Cavity  at  Different  L/R 
Ratios. 


as  well  as  the  reflected  radiation  would  be  attenuated.  A  separate  linearity  check  was  made 
with  each  experimental  set-up,  and  the  data  are  shown  in  Tables  1  and  2,  (pages  53  and  54). 

2.4b  First  Reflection  Cavity  Tests 

The  procedures  used  to  calibrate  the  integrating  sphere  reflectometer  were  based  on  pre- 
vious studies  by  Kelly  and  Moore  [  2  ].  They  verified  experimentally  the  Gouf f  e  [  3 1  equation 
relating  the  reflectance  of  a  shallow  cavity  to  its  geometry  and  reflectance  of  its  walls  by 
measuring  the  reflectance  of  shallow  cavaties  of  known  properties  with  a  General  Electric 
recording  spectrophotometer. 

If  a  cavity  is  used  that  has  a  diffusely  reflecting  bottom  and  completely  absorbing 
walls,  only  the  flux  that  escapes  through  the  opening  on  the  first  reflection  is  available 
for  measurement.  Such  a  cavity  is  shown  in  Figure  5a.  The  depth  (L)  of  the  cavity  is  adjust- 
able, so  that  L/R  ratios  (R  is  the  radius  of  the  cavity  opening)  can  be  varied  from  0.0  to 
2.0.  The  reflectance  of  the  cavity  will  then  be 

Pc  =  Pwf  (3) 

where  Pc  is  the  reflectance  of  the  cavity,  Pw  is  the  diffuse  reflectance  of  the  cavity  bottom, 
and  f  is  the  diffuse  configuration  factor  between  the  cavity  bottom  and  opening.  For  a  small 
spot  at  the  center  of  the  cavity  bottom, 


1  +  (L/R)2 


(4) 


The  tests  with  the  first  reflection  cavity  are  particularly  appropriate  for  calibration 
of  the  integrating  sphere  reflectometer  for  two  reasons:   (l)  Equation  (4)  has  been  experi- 
mentally verified  with  an  instrument  of  proven  reliability,  and  (2)  the  geometric  distribution 
of  reflected  flux  varies  from  a  solid  angle  of  2 rr steradians  at  L/R  =  0  (cone  with  central 
plane  angle  of  180°)  to  about  0.84  steradian  at  L/R  =  2  (cone  with  central  plane  angle  of  60°) 
as  illustrated  in  Figure  6.  Failure  of  an  integrating  sphere  to  give  measured  F  values  that 
agree  with  the  f  values  computed  from  equation  (4)  could  result  from  either  deviation  of  the 
reflectance  of  the  cavity  bottom  from  perfect  diffuseness,  or  improper  viewing  conditions  for 
the  detector. 

Figure  6  shows  the  areas  of  the  sphere  that  are  directly  illuminated  by  the  first  reflec- 
tion cavity  at  different  L/R  ratios.  It  is  apparent  that  at  low  L/R  ratios  the  detector  port  . 
is  directly  illuminated.  If  the  detector  views  only  a  portion  of  the  sphere  wall,  the  frac- 
tion of  the  viewed  area  that  is  directly  illuminated  will  vary  with  L/R,  particularly  at  high 
values  of  L/R.   Consequently  the  detector  response  will  depend  not  only  on  p^f  but  also  on 
the  area  of  the  sphere  directly  illuminated  by  the  cavity. 

The  area  viewed  by  the  detector  can  be  varied  by  use  of  different  accessories  between 
the  detector  and  the  detector  port.  Figure  7  is  a  schematic  diagram  of  the  detector  assembly, 
which  consists  of  four  parts;  (l)  an  opal  glass  viewing  window,  (2)  sleeving,  which  may  be 
either  white  or  black,  (3)  a  light  pipe,  and  (4)  a  photomultiplier  and  cover.  The  space  be- 
tween the  detector  port  and  detector  can  be  occupied  by  any  one  or  more  of  the  elements 
1,  2,  and  3  above.  As  an  alternative,  a  series  of  baffles  may  be  introduced  into  the  passage 
between  the  detector  and  port  to  reduce  the  area  viewed. 

Measurements  were  made  of  the  reflectance  of  the  first  reflection  cavity  with  a  paper 
of  3M  Velvet  yellow  applied  to  the  base,  and  with  L/R  settings  of  0.0,  0.50,  0.75,  1.00,  1.50 
and  2.00.   Several  different  detector  viewing  configurations,  as  illustrated  in  Figure  8,  were 
used.  The  measured  reflectance  at  each  value  of  L/R  was  divided  by  the  reflectance  at  L/R=0 
to  obtain  the  F  values  shown  in  Table  1.  The  difference  between  the  computed  values,  f,  and 
measured  values,  F,  is  plotted  as  a  function  of  L/R  in  Figures  9,  10,  and  11. 

The  data  show  that  the  configurations  with  the  opal  glass  viewing  cap  (trials  105  and 
106),  which  enable  the  detector  to  view  the  entire  sphere  wall,  give  high  values  of  F.  The 
configurations  that  restrict  the  field  of  view  with  collimators  to  a  small  area  opposite  t] 
detector  (trials  108,  109,  110),  and  with  black  sleeves  only  (trials  104  and  107),  give  low 
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Figure  8.  Detector  Con- 
figurations Used  in 
First  Reflection  Cavity 
Tests . 


Figure  9.  Difference  Between  Measured  F 
and  Computed  f ,  Plotted  as  a  Function 
of  L/R,  for  Detector  Configurations 
With  Black  Sleeving,  101-Black  Sleeve 
with  Light  Pipe;  104-Black  Sleeve  Only; 
107-Long  Black  Sleeve,  0.075  Inch  Beam. 


L/R 


Figure  10.  Difference  Between  Measured 
F  and  Computed  f ,  Plotted  as  a  Function 
of  L/R,  for  Detector  Configurations  With 
Collimators,  108  -  5-Inch  Collimator, 
0.075  Inch  Beam;  109  -  ?-Inch  Collimator, 
0.075  Inch  Beam;  110  -  Same  as  109, 
Except  EME  5685  Photomultiplier  Instead 
of  IP28. 
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Figure  11.  Difference  Between  Measured 
F  and  Computed  f ,  Plotted  as  a  Function 
of  L/R,  for  Detector  Configurations 
with  White  Sleeving,  105  -  White  Sleeve, 
Opal  Glass  View  Cap  and  Light  Pipe; 
106  -  Same  as  105,  but  0.075  Inch  Beam; 
111  -  White  Sleeve,  Light  Pipe,  0.075 
Inch  Beam. 


values  of  F.  Best  agreement  was  obtained  in  trials  101  (black  sleeve,  light  pipe)  and  111 
(white  sleeve,  light  pipe,  0.075  in.  beam). 

The  results  of  these  tests  suggested  that  the  integrating  sphere  was  only  marginally 
suitable  for  the  Intended  use.  It  has  not  been  established  whether  this  is  due  to  (l)  the 
relatively  low  reflectance  of  the  sphere  coating,  about  0.82  as  applied,  (2)  variations  in 
reflectance  over  the  surface  of  the  sphere,  (3)  the  location  of  the  detector  port,  or  (4) 
the  lack  of  baffling  between  the  detector  port  and  specimen. 

2.4c  Lined  Cavity  Tests 

If  a  cavity  is  used  which  has  diffusely  reflecting  bottom  and  walls,  as  illustrated  in 
Figure  5b,  a  smaller  fraction  of  the  incident  flux  will  be  absorbed,  and  the  reflected  flux 
will  be  more  nearly  diffuse.  The  experimental  data  of  Kelly  and  Jfoore  [2]  agreed  well  with 
those  computed  from  the  Gouff €  equation  for  the  reflectance  of  the  cavity,  pc  ,  in  terms  of 
the  L/R  ratio  and  the  reflectance  of  the  cavity  walls,  pN . 


(l-pH)  [l  +  pH(f  -f)]         f-pH(f-f) 

Pc  ~  1  ~  -   Pw 


(5) 


1-Ph(1  -!)  1  -  Ph(1  -A) 


a      1 1        a 

where  7  =  ,..    ,   r/n\     and  f  =  1  +  (l/r)2  '  A  is  tne  ratio  of  tne  opening  of  the  cavity  to 

the  total  wall  area  of  the  cavity,  opening  included,  and  f  is  the  fraction  of  the  incident 
flux  leaving  the  cavity  on  the  first  reflection  (equation  4)»  Setting  H  =  pc/pM  gives 

H  -  f 
Pn  = (6) 

(H  -  f )  +  f   (1  -  H) 

Equations  (5)  and  (6),  like  equation  (4),  are  based  on  center-point  illumination  and  on 
the  assumption  that  the  walls  and  bottom  of  the  cavity  are  perfectly  diffuse  reflectors.  The 
laser  beam  gives  a  good  approximation  of  the  center  point  illumination,  as  the  diameter  of 
the  beam  was  only  about  1.9mm,  while  the  diameter  of  the  cavity  was  32mm. 

If  the  values  measured  with  the  first  reflection  cavity  agree  with  the  computed  f  values, 
the  lined  cavity  values  should  yield  accurate  absolute  reflectances.  The  effect  of  uncer- 
tainties in  L/R  and  the  measured  reflectance  ratio,  Qc/0^,  on  the  accuracy  of  the  computed 
absolute  reflectance,  pH ,  is  shown  in  Figures  12  and  13 •  The  data  indicate  that  the  method 
should  give  reasonable  precision,  even  for  reflectances  as  low  as  0.50.  The  method  is  con- 
venient, since  a  cavity  is  easily  prepared.  If  the  reflectometer  gives  accurate  values  of 
F  with  the  first  reflection  cavity  and  pHwith  the  lined  cavity,  it  can  be  assumed  that 
accurate  values  will  be  obtained  with  diffusely  reflecting  specimens  of  almost  any  geometric 
distribution  of  reflected  flux,  provided  there  is  negligible  specular  reflection  back  out 
the  entrance  port. 

The  measured  values  of  F  obtained  in  trials  101,  104,  and  111,  reported  above,  gave 
the  best  agreement  with  computed  f  values .  Hence  evaluation,  of  pM ,  the  absolute  reflectance 
of  the  3M  Velvet  yellow  and  green  papers,  was  attempted  by  the  shallow  cavity  approach,  using 
these  same  detector  configurations.  The  data  for  the  yellow  paper  are  reported  in  Table  2 
as  trials  102,  103,  and  112 ,  with  the  viewing  configurations  reported  in  trials  101,  104, 
and  112,  respectively,  and  for  the  green  paper  as  trials  301  for  the  first  reflection  cavity 
and  302  for  the  lined  cavity,  using  the  same  viewing  configuration  as  in  trial  111.  The 
data  were  reduced  by  use  of  the  following  equation. 


P„  =  H~F  (7) 
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Figure  13.  Effect  of  Uncertainty  in 
Measured  Reflectance  Ratio  Qc/Qw,  on 
Accuracy  of  pw  Evaluated  by  Shallow 
Cavity  Method. 
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Figure  14-  Absolute  Reflectance  of 
3M  Velvet  Yellow  Paper  Computed  from 
Data  of  Trails  101  and  102,  Top,  and 
103  and  104,  Bottom,  Plotted  as 
Functions  of  L/R  Ratio. 
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Figure  15.  Absolute  Reflectance  of 
3M  Velvet  Yellow  Paper,  Computed 
from  Data  of  Trials  111  and  112, 
Top,  and  of  3M  Velvet  Green,  Com- 
puted from  Data  of  Trials  301  and 
302,  Bottom,  Plotted  as  Functions 
of  L/R  Ratio. 
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Figure  16.  Experimentally  Determined  F 
Values  and  Absolute  Reflectance  pH ,  of 
3M  Velvet  Yellow  Paper,  Computed  by- 
Use  of  Equation  (7),  Trials  101-102. 
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Figure  17-  Experimental 3y  Determined  F 
Values  and  Absolute  Reflectance,  pH  ,  of 
3M  Velvet  Yellow  Paper,  Computed  by 
Use  of  Equation  (7),  Trials  103-104- 
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Figure  18.  Experimentally  Determined  ] 
Values  and  Absolute  Reflectance,  pw , 
3M  Velvet  Yellow  Paper,  Computed  by 
Use  of  Equation  (7),  Trials  111-112. 
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Figure  19-  Experimentally  Determined  F 
Values  and  Absolute  Reflectance  of  3M 
Velvet  Green  Paper,  Computed  by  Use 
of  Equation  (7),  Trials  301-302. 
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[Equation  (7)  is  a  modification  of  equation  (6)  in  which  a  distinction  is  made  between  F  and 
If.  F  is  used  in  the  equation  where  flux  from  the  first  reflection  only  is  considered,  and 
jf  where  the  flux  from  multiple  reflections  is  involved.  Ideally,  F  =  f  for  a  completely 
.diffuse  cavity  wall  and  an  ideal  integrating  sphere  reflectometer.  However,  the  measured  F 
'deviated  somewhat  from  the  computed  f  in  nearly  every  case.  H,  F,  f  and  a/A  are.  as  previ- 
jously  defined. 

The  data  were  reduced  by  three  procedures: 

(1)  Computed  values  of  f  were  used  for  both  f  and  F. 

(2)  Measured  values  of  F  were  used  for  both  f  and  F. 

(3)  Measured  values  of  F  and  computed  values  of  f  were  used  where  applicable. 

The  results  of  the  computations  are  given  in  Table  3  for  the  original  data  obtained  in 
! trials  101-102,  103-104,  111-112  and  301-302,  and  are  plotted  as  functions  of  L/R  ratio  in 
Figures  14  and  15.  The  results  obtained  by  procedure  3,  together  with  the  corresponding 
lvalues  of  f  and  F,  are  plotted  as  functions  of  L/R  ratio  in  Figures  16,  17,  18  and  19. 

Procedure  3  appears  to  give  the  most  consistent  results.  The  over-all  scatter  is  least, 
jand  the  three  curves  for  p  of  the  3M  Velvet  yellow  paper  computed  by  this  procedure  nearly 
[coincide.  The  reflectance  of  the  3M  Velvet  paper,  relative  to  MgO,  was  O.85O  at  0.632u  as 
measured  with  the  General  Electric  recording  spectrophotometer.  Converting  this  reflectance 
|to  absolute  on  the  basis  of  the  absolute  reflectance  of  MgO  of  0.985  at  0.632u  as  reported 
iby  Edwards  et  al  [4],  we  obtain  a  value  of  0.837.  The  corresponding  value  for  the  reflect- 
jance  of  the  3M  Velvet  green  paper,  relative  to  MgO,  is  O.385  at  0.632u,  which  converts  to 
'0.379  on  an  absolute  scale. 

All  of  the  absolute  reflectances  computed  from  data  taken  at  L/R  values  of  0.5  appear 
ito  be  in  error.  The  values  are  high  as  computed  by  methods  1  and  2,  and  low  as  computed  by 
method  3.  There  is  also  a  tendency  for  the  absolute  reflectances  computed  from  data  taken 
I at  L/R  values  of  2.0  to  be  somewhat  low.  The  data  computed  by  procedure  3  appear  to  be 
Ithe  most  consistent.  The  12  values  obtained  by  this  procedure  at  L/R  values  of  0.75,  1.0, 
|1. 5  and  2.0,  average  0.842  with  a  standard  deviation  of  0.006,  which  is  an  excellent  check 
I  of  the  value  of  0.838  computed  from  the  reflectance  measured  with  the  General  Electric  re- 
cording spectrophotometer.  The  corresponding  average  value  for  the  absolute  reflectance 
lof  the  3M  Velvet  green  paper  at  0.632u,  however,  is  0.339,  which  is  significantly  lower  than 
the  value  of  O.38O  computed  from  the  reflectance  measured  with  the  General  Electric  record- 
ing spectrophotometer. 

The  flux  from  the  first  reflection,  ioPw^,  is  a  major  portion  of  the  total  flux  reflect- 
jed  from  a  shallow  cavity.  K,  the  fraction  of  the  total  flux  reflected  on  the  first  reflec- 
tion, can  be  computed  for  a  perfect  diffuser  as 

K  =  ^  = 


p  f 

Kw 

*-P„*U-t> 

>c 

*  -«.(}-*) 

The  fraction  K  is  plotted  as  a  function  of  wall  reflectance,  p  ,  for  different  L/R 
ratios  in  Figure  20.  The  data  show  that  for  low  p  ,  K  becomes  very  large,  that  K  decreases 
with  increasing  L/R  for  any  p  ,  and  that  K  decreases  markedly  with  an  increase  in  L/R  at 
high  p  .  As  a  first  approximation,  the  fraction  K  of  the  flux  reflected  from  a  shallow 
I  cavityVLll  be  distributed  uniformly  over  the  solid  angle  subtended  by  the  cavity  opening 
'from  the  center  of  the  incident  beam,  and  the  fraction  (1  -  K)  of  the  flux  will  be  dis- 
tributed uniformly  over  the  entire  hemisphere.  Hence  the  geometric  distribution  of  the  flux 
reflected  from  a  shallow  lined  cavity  will  be  intermediate  between  that  from  a  first  reflec- 
tion cavity  and  that  of  a  perfect  diffuser. 

No  detailed  explanation  of  the  shape  of  the  curves  in  Figures  14  through  19  will  be 
attempted  at  this  time.  The  data  are  preliminary,  and  are  subject  to  errors  that  will  be 
evaluated  as  the  work  progresses. 

2.4d  Appended  Sphere  Tests 

The  absolute  reflectance  of  a  coating  material  can  also  be  evaluated  by  the  appended 
sphere  meth.0d.3_/  The  diffusely  reflecting  coating  to  be  evaluated  is  applied  to  the  inside 


2/  Much  of  the  theory  for  the  appended  sphere  method  was  supplied  by  David  D.  Goebel 
of  the  Photometry  and  Colorimetry  Section  of  the  National  Bureau  of  Standards. 
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Figure  20.  Ratio  of  First  Reflected  Flux 
to  Total  Flux,  K,  Reflected  by  a  Shallow 
Cavity  for  Different  L/R  Ratios,  as  a 
Function  of  Wall  Reflectance,  pw. 
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Figure  21.  Relation  of  Reflectance,  p,  to 
Measured  Reflectance  Ratio,  Qc/Qf,  in  the 
Appended  Sphere  Method,  for  Different 
Sphere  Geometries,  d/D. 
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ance by  the  Appended  Sphere  Method, 
Plotted  as  a  Function  of  Reflectance, 
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Figure  23.  Effect  of  Change  in  Sphere  Geometry, 
d/D  on  Reflectance  by  the  Appended  Sphere 
Method,  Plotted  as  a  Function  of  Reflectance, 
p,  for  Different  Sphere  Geometries. 
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surface  of  a  small  sphere  having  a  single  hole  of  the  same  diameter  as  the  specimen  port  on 
the  integrating  sphere,  and  to  a  flat  specimen  of  any  convenient  size.  The  reflectance  of 
(the  sphere,  Q„,  and  of  the  flat  specimen,  Q„,  are  measured.  These  reflectances  can  be 
measured  relative  to  any  convenient  standard  or  to  each  other. 

The  absolute  reflectance  of  the  coating,  p,  is  computed  from  the  measured  reflectance 
I  ratio,  Qs/Qf,  by  use  of  the  following  equation: 


1  - 


Qo/Qf 

P  -  —  (9) 

1-q 

where  q  is  the  fraction  of  incident  flux  that  leaves  the  sphere  on  the  first  reflection 
computed  as: 

q 


-Uw1-^]  <10> 


where  d  is  the  diameter  of  the  opening  in  the  appended  sphere,  and  D  is  the  diameter  of  the 
appended  sphere. 

The  relation  of  p  to  Qo/Qf.  for  different  values  of  d/D  is  shown  in  Figure  21.  These 
curves  indicate  that  as  the  reflectance  p  of  the  coating  decreases,  the  reflectance  ratio 
Qo/Qf  decreases  very  rapidly,  particularly  for  smaller  d/D.   The  conclusion  is  that  the 
method  is  not  very  precise  for  measuring  absolute  reflectance  of  low-reflectance  coatings. 
However,  the  appended  sphere  method  does  have  inherently  high  accuracy  for  evaluating  the 
absolute  reflectance  of  highly  reflecting  diffuse  coatings,  as  is  shown  in  Figures  22  and 
23.  Figure  22  shows  the  uncertainty  in  p  produced  by  an  uncertainty  in  the  measured  reflect- 
ance ratio,  Qc/Qf,  for  different  d/D  ratios  as  a  function  of  p.  Figure  23  shows  the  cor- 
responding uncertainties  in  p  due  to  uncertainties  in  d/D.  As  an  example,  for  d/D  of  0.250, 
if  Qo/Qf.  is  uncertain  by  1  percent  and  d/D  is  also  uncertain  by  1  percent,  the  uncertainty 
in  p  at  0.98  is  only  0.2  percent.  However,  if  p  is  only  0.80,  the  corresponding  uncertainty 
in  p  increases  to  about  4o4  percent. 

The  appended  sphere  method  differs  from  the  shallow  cylindrical  cavity  method  in  two 
important  respects:   (l)  there  is  no  theoretical  verification  for  the  appended  sphere 
method,  and  (2)  the  flux  reflected  from  the  appended  sphere  has  a  geometric  distribution 
that  more  closely  approximates  that  from  a  perfect  diffuser. 

Several  spheres  were  prepared  for  use  with  the  appended  sphere  method,  but  no  measure- 
ments were  completed  during  the  report  period. 

2.5  Relative  Reflectance  Measurements 

The  laser— source  integrating  sphere  reflectometer  should  be  suitable  for  making  rela- 
tive reflectance  measurements,  even  if  errors  are  present  in  the  measured  absolute  reflect- 
ance. The  3M  Velvet  yellow  paper  was  selected  as  a  comparison  standard,  and  the  relative 
reflectance  of  a  number  of  materials  was  measured.   The  data  are  given  in  Table  4;  and  com- 
pared to  similar  values  obtained  on  the  same  specimens  with  the  General  Electric  recording 
spectrophotometer  or  the  Cary  Model  14  spectrometer,  both  at  0.632|i.   The  data  are  pre- 
liminary, and  no  detailed  analysis  will  be  attempted  at  this  time. 

2.6  Modifications  to  Reflectometer 

Because  of  the  somewhat  disappointing  results  obtained  with  the  integrating  sphere  re- 
flectometer as  originally  designed  and  with  the  3M  Velvet  white  coating,  the  sphere  was 
modified  by  adding  a  second  detector  port,  135°  from  the  specimen  port  in  the  plane  of  the 
entrance  port.  Because  of  the  location  of  the  flange  connecting  the  two  halves  of  the 
sphere,  it  is  also  possible  to  rotate  the  lower  half  of  the  sphere  so  that  the  original 
|  detector  port  is  in  a  plane  normal  to  the  plane  of  the  entrance  port.  The  3M  coating  was 
replaced  by  BaSO,  coating,  having  appreciably  higher  reflectance.  It  is  expected  that  the 
new  sphere  coating  and  detector  port  will  appreciably  reduce  the  errors  that  were  present 
in  the  original  design. 
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3.  Ellipsoidal  Mirror  Reflectometer 

3.1  Background  -  Need  for  Instrument 

Normal  spectral  emittance  measurements  of  specimens  at  and  near  room  temperature  are  not 
easy  to  make.  The  total  amount  of  energy  emitted  at  such  temperatures,  even  by  a  blackbody 
radiator,  is  small,  and  when  further  attenuated  by  a  monochromator,  the  amount  is  so  small 
that  special  techniques  are  required  to  measure  it.  On  the  other  hand,  energy  considerations 
do  not  seriously  affect  reflectance  measurements.  As  long  as  the  temperature  of  the  specimen 
remains  approximately  constant,  the  fraction  of  the  incident  flux  that  is  reflected  is  inde- 
pendent of  the  incident  flux  density.  Hence  we  can  measure  spectral  reflectance  and  compute 
spectral  emittance  by  use  of  Kirchhoff 's  law. 

There  is  a  more  urgent  and  specific  need  for  spectral  reflectance  data  in  the  current 
program.  The  laser-source  reflectometer  described  in  Section  2  of  this  report  measures  re- 
flectance only  at  those  few  wavelengths  at  which  continuous  wave  lasers  operate,  but  provides 
some  spectral  reflectance  data  at  high  temperatures.   It  is  fundamentally  incapable  of  provid- 
ing the  continuous  spectral  reflectance  curve  needed  for  heat-transfer  computations. 

Available  data  indicate  that,  while  spectral  emittance  or  reflectance  of  a  material  at  any: 
one  wavelength  may  change  significantly  with  temperature,  the  general  shape  of  the  spectral 
reflectance  or  emittance  curve  .does  not,  provided  the  material  does  not  undergo  any  phase 
change  or  irreversible  change  in  physical  properties  as  a  result  of  heat  treatment.  Hence  it 
is  possible  to  construct  a  spectral  reflectance  or  emittance  curve  for  material  from  accurate 
data  obtained  at  a  few  wavelengths  with  the  laser  source  reflectometer,  by  interpolating 
values  at  intermediate  wavelengths  on  the  basis  of  the  shape  of  the  room  temperature  curve. 
For  the  reflectance  to  be  useful  for  computing  normal  spectral  emittance  by  Kirchhoff 's  law, 
or  for  interpolation  of  laser-source  reflectometer  data,  it  must  be  measured  under  equivalent  - 
optical  conditions  —  normal  illumination  and  hemispherical  viewing,  or  the  optical  equivalent 
diffuse  illumination  and  normal  viewing. 

A  number  of  different  instruments  have  been  designed  to  measure  reflectance  under  condi- 
tions approximating  normal  illumination  and  hemispherical  viewing.  These  include  the  Coblentz. 
hemisphere,  which  approximates  normal  illumination  and  hemispherical  viewing,  the  integrating 
sphere  reflectometer,  which  can  be  used  to  approximate  either  of  the  conditions  mentioned 
above,  and  the  hohlraum  reflectometer  and  Janssen's  modification  of  the  Coblentz  hemisphere 
reflectometer,  which  approximate  diffuse  illumination  and  hemispherical  viewing.  Each  of  thes 
has  advantages  and  disadvantages,  but  none  was  deemed  entirely  suitable  for  our  measurements. 

In  an  attempt  to  design  an  instrument  that  would  provide  data  of  the  desired  reliability,  • 
at  wavelengths  out  to  at  least  15  microns,  the  ellipsoidal  mirror  reflectometer  was  conceived. 
Basically,  the  instrument  is  a  modification  of  the  Coblentz  hemisphere,  in  which  an  ellipsoida 
mirror  replaces  the  hemispherical  mirror.  A  narrow  beam  of  radiation  from  the  monochromator 
is  focused  through  a  small  hole  in  the  ellipsoidal  mirror  onto  the  specimen,  which  is  located 
at  the  first  focus  of  the  ellipsoid.  The  ellipsoidal  mirror  focuses  the  radiation  reflected 
by  the  specimen  onto  the  detector,  which  is  located  at  the  second  focus  of  the  ellipsoid,  some 
17  inches  below  the  specimen. 

Advantages  of  the  ellipsoidal  mirror  reflectometer  include: 

A.  The  detector  views  a  solid  angle  having  a  central  plane  angle  of  only  about  48°  rather 
than  a  hemisphere,  as  in  the  Coblentz  method,  hence  errors  due  to  variations  in  angular  sensi- 
tivity of  the  detector  are  reduced. 

B.  The  ellipsoidal  mirror  produces  an  image  that  is  much  freer  from  optical  aberrations 
than  is  the  case  with  a  hemispherical  reflector. 

C.  The  separation  of  sample  and  detector  by  some  17  inches  permits  the  specimen  to  be 
heated  or  cooled  without  adversely  affecting  the  detector. 

D.  Diaphragms  at  the  first  focal  plane  can  be  used  to  selectively  block  out  or  pass  ra- 
diation reflected  into  known  solid  angles,  hence  measurements  can  be  made  of  the  geometric 
distribution  of  reflected  radiation. 
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E.  The  reflectance  of  the  mirror  coating  remains  high  throughout  the  infrared.  Most 
ntegrating  sphere  coatings  have  low  reflectance  beyond  about  2. 5  n. 


F.  The  instrument  can  be  used  for  absolute  as  well  as  comparative  reflectance  measure- 
ents. 

The  principal  disadvantages  of  the  ellipsoidal  mirror  reflectometer  are: 

A.  The  specimen  and  its  support  block  a  small  portion  of  the  optical  path  at  the  first 
focal  plane. 

B.  The  mirror  has  a  linear  magnification  factor  of  about  5,  hence  a  large-area  detector 
Ls  required,  and  the  illuminated  area  on  the  specimen  must  be  kept  small. 

3.2  Previous  Work 

,2a  Design  of  Equipment 

Equipment  for  measuring  spectral  reflectance  under  conditions  approximating  normal  illu- 
mination and  hemispherical  viewing,  over  the  wavelength  range  of  approximately  1  to  15  microns, 
was   designed  and  built,  as  described  in  WADC  TR  59-510  pts.   II  [  52]and  III T  53  ]. 

The  reflectometer  consisted  of  six  distinct  parts:   (1)  a  source  assembly,  (2)  a  monochro- 
imator,  (3)  transfer  optics,  (4)  an  ellipsoidal  mirror,  (5)  a  detector,  and  (6)  an  amplifier 
(and  recorder.  The  source  assembly,  (l),  consisted  of  a  silicon-carbide  element,  a  chopper  and 
la  pair  of  mirrors,  one  plane  and  one  spherical,  which  focused  the  chopped  flux  from  the  hot 
element  onto  the  entrance  slit  of  the  monochromator  (2).  The  monochromator,  (2),  which  was 
equipped  with  a  sodium  chloride  prism,  dispersed  the  flux.  The  transfer  optics,  i3) ,   directed 
and  focused  the  chopped  and  dispersed  flux,  through  a  small  hole  in  the  ellipsoidal  mirror, 
onto  the  specimen,  which  was  located  at  the  first  focus  of  the  ellipsoid.  The  ellipsoidal 
mirror  (4)  collected  the  flux  reflected  by  a  specimen  located  at  its  first  focal  point  and 
focused  it  on  the  detector  located  at  the  second  focal  point.  A  Golay-cell  detector  (5)  was 
mounted  on  a  device  which  permitted  it  to  be  positioned  at  either  the  first  focal  point  or 
second  focal  point  of  the  mirror.  The  amplifier  and  recorder  (6),  made  a  permanent  record  of 
the  amplified  signal  from  the  detector.   The  reflectometer  was  designed  to  measure  absolute 
reflectance  by  the  following  procedure. 

The  incident  flux  was  measured  with  the  detector  at  the  first  focus  of  the  ellipsoid;  the 
detector  was  then  moved  to  the  second  focus  of  the  ellipsoid  and  the  specimen  was  mounted  at 
the  first  focus,  and  the  reflected  flux  was  measured.  The  ratio  of  the  reflected  flux  to  in- 
cident flux,  corrected  for  mirror  losses,  was  taken  as  the  absolute  spectral  reflectance  of 
the  specimen  at  the  wavelength  of  the  incident  flux. 

3.2b  Golay-cell  Detector  Problems 

A  large-area  Golay-cell  detector  was  used  in  the  early  work.   Its  sensitive  area  was  one 
centimeter  in  diameter.  When  used  with  the  sensitive  diaphragm  in  a  horizontal  position  the 
detector  proved  to  be  highly  microphonic.  Elaborate  precautions  were  taken  to  isolate  the 
detector  from  structure-borne  noise  by  mounting  it  on  a  seismic  table  which  was  separated  from 
the  rest  of  the  equipment,  and  to  isolate  it  from  air-borne  noise  by  enclosing  the  entire 
optical  path  from  the  exit  port  of  the  monochromator  to  the  detector  in  acoustical  insulation. 
These  precautions  reduced  the  noise  to  a  point  where  useful  measurements  could  be  made. 

Measurements  with  this  early  design  of  the  reflectance  of  both  flowers  of  sulfur  and  of 
a  rhodium  mirror  gave  values  of  reflectance  as  low  as  0.50.  The  reflectance  of  these  materi- 
als at  the  wavelengths  at  which  measurements  were  made  was  known  to  be  greater  than  0.90. 
Mirror  losses  were  computed  to  be  appreciably  less  than  0.05  in  the  case  of  sulfur  and  0.03 
in  the  case  of  the  rhodium  mirror.  Obviously  the  measured  reflectances  were  grossly  in  error. 

One  of  the  objectives  of  the  work  with  the  instrument  was  to  measure  the  reflectance  of 
specimens  at  temperatures  above  and  below  room  temperature.  A  specimen  heater,  capable  of 
operation  up  to  800°K,  was  designed  and  built.  When  a  test  was  made  of  a  heated  specimen,  the 
Golay-cell  detector  failed  because  of  the  background  radiation  when  the  specimen  reached  a 
temperature  of  about  350°  K. 
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Because  of  these  difficulties  encountered  with  the  Golay-cell  detector,  the  decision 
was  reached  that  a  different  type  of  detector  was  needed. 

3.3  Current  Work 

3.3a  Detector  Developments 

3.3a-l  Thermopile  Detectors 

Two  custom-built  thermopile  detectors  with  a  viewing  area  of  one  square  centimeter  were 
obtained  and  tested.  Tests  of  signal-to-noise  ratio,  sensitivity  and  linearity  in  the  1  to 
15  micron  wavelength  range  indicated  that  they  should  be  satisfactory  for  the  intended  purpose 

A  preliminary  measurement  of  the  reflectance  of  flowers  of  sulfur  at  5m-  gave  an  uncor- 
rected reflectance  of  about  0.50,  which  is  obviously  in  error.   It  was  thought  that  some  of 
the  losses  might  be  caused  by  aberrations  of  the  ellipsoidal  mirror,  which  resulted  in  some 
of  the  flux  reflected  from  the  specimen  being  focused  outside  the  sensitive  area  of  the 
detector.   Photographs  were  taken  of  (l)  the  beam  incident  on  the  specimen,  and  (2)  the  beam 
incident  on  the  detector  after  reflection  from  the  specimen.   By  careful  adjustment  of  the 
position  of  the  specimen  and  detector  relative  to  the  ellipsoidal  mirror,  it  was  possible 
to  obtain  an  image  at  the  detector  of  the  illuminated  spot  on  the  specimen  that  fell  entire- 
ly within  the  sensitive  area  on  the  detector.  However,  this  careful  adjustment  did  not 
appreciably  increase  the  measured  (uncorrected)  reflectance  of  flowers  of  sulfur  at  5u. 

3.3a-2  Spatial  Sensitivity  of  Detectors 

Further  investigation  revealed  that  both  the  Golay  and  thermopile  detectors  gave  varied 
response  over  their  sensitive  areas.   Quantitative  tests  were  made  to  show  the  magnitude  of 
this  effect. 

The  Golay-cell  detector  was  mounted,  with  its  sensitive  area  in  a  vertical  plane,  on  an 
automatically-driven  micrometer  head  that  could  move  it  horizontally  in  the  plane  of  the 
sensitive  area  at  a  rate  of  about  0.08"  per  minute.  A  stationary  aperture  stop,  having  a 
circular  opening  l/l6"  in  diameter,  was  mounted  directly  in  front  of  the  detector.  The  de- 
tector was  illuminated,  through  the  aperture  stop,  by  filtered  radiation  from  the  Globar 
source,  chopped  at  10  cycles  per  second.  Three  different  filter  cc.      ;ns  were  used,  to 
give  two  relatively  narrow  bands  centered  at  about  2  and  9  microns  respectively,  and  one 
covering  the  range  of  15  to  20  microns. 

Two  scans  across  the  sensitive  area  of  the  cell  were  made  with  radiation  passing  each 
of  the  filter  combinations,  along  two  diameters  of  the  cell  normal  to  each  other,  designated 
A-A'  and  B-B'  respectively.  Results  of  the  scans  are  shown  in  Figures  24  and  25.  Each 
curve  represents  the  relative  response  of  the  detector-amplifier  combination  to  the  radiant 
energy  passing  through  the  aperture,  as  a  function  of  the  distance  along  the  scan.  Differ- 
ent gain  settings  were  used  for  the  different  spectral  regions,  hence  the  different  curves 
are  not  directly  comparable. 

The  curves  in  Figure  24  show  that  along  the  diameter  A-A1  an  area  of  highest  sensitivity 
exists  2-3  mm  from  the  right  side  of  the  sensitive  area,  and  a  second,  much  lower  peak  in 
sensitivity  occurs  near  the  center.  These  results  were  not  expected,  and  there  appears  to  be 
nothing  in  the  basic  design  of  the  detector  to  account  for  them.  Apparently  the  absorptance 
of  the  absorbing  membrane  varies  over  its  surface.   The  curves  in  Figure  25  show  that  along 
the  diameter  B-B'  the  sensitivity  is  symmetrical,  with  two  peaks  each  about  3mm  from  tJ 
edge  of  the  sensitive  area,  with  slightly  less  sensitivity  at  the  center.  There  is  a  short- 
range  periodicity,  with  a  minor  peak  about  every  millimeter.  The  curves  in  Figure  25  are 
more  nearly  what  would  be  expected  from  the  basic  design  of  the  cell. 

Similar  measurements  were  made  of  the  variation  in  sensitivity  over  the  sensitive  area 
of  the  new  thermopile  detector  that  was  referred  to  above.  In  this  case  the  Globar  source 
was  replaced  by  a  tungsten  filament  lamp,  and  no  filters  were  used.  A  diffusing  screen  was 
placed  immediately  in  front  of  the  detector  for  some  of  the  scans.   The  diffusing  screen  was 
a  sodium  chloride  window  that  had  been  depolished  on  the  side  away  from  the  detector.  Dif- 
fusing screen  no.  2  was  5  mm  thick,  and  one  surface  had  been  ground  with  9.5n  abrasive. 
Screen  no.  3  was  2.5mm  thick,  and  one  surface  had  been  ground  with  50u  abrasive. 

The  thermopile  consists  of  ten  plates,  each  approximately  2  x  5mm  in  size,  arranged  in 
two  columns  of  five  rows  each  to  form  a  sensitive  area  one  centimeter  square.   A  thermocouple 
was  attached  to  the  back  of  each  plate,  and  the  10  thermocouples  were  connected  in  series 
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Figure  24.  Spatial  Sensitivity  of 
the  Golay-Cell  Detector  in  the 
A-A1  Direction.   The  curves  show 
the  variation  in  response  to  a 
beam  of  constant  flux,  1.6  mm 
(l/l6  inch)  in  diameter,  as  the 
beam  is  scanned  across  the  de- 
tector. Flux  of  three  different 
wavelengths  was  used,  as  indicated. 


DIFFUSING  SCREEN  NO.  2 
B'  —  B 

Figure  26.  Spatial  Sensitivity  of  the  Ther- 
mopile Detector  in  the  A-A*  and  B-B1  Di- 
rections to  Chopped  Incident  Flux.   The 
curves  show  the  variation  in  response  to 
a  beam  1.6  mm  (l/l6  inch)  in  diameter  of 
chopped  incident  flux  from  a  constant 
tungsten  source,  as  the  beam  is  scanned 
across  the  detector.  The  effect  of  two 
different  diffusing  screens  in  reducing 
spatial  sensitivity  can  be  seen. 


Figure  25.  Spatial  Sensitivity  of 
the  Golay-Cell  Detector  in  the 
B-B'  Direction.   The  curves  show 
the  variation  in  response  to  a 
beam  of  constant  flux,  1.6  mm 
(l/l6  inch)  in  diameter,  as  the 
beam  is  scanned  across  the  de- 
tector. Flux  of  three  different 
wavelengths  was  used,  as  indicated. 


Figure  27.   Spatial  Sensitivity  of  the  Ther- 
mopile Detector  in  the  A-A'  and  B-B'  Di- 
rections to  Unchopped  Incident  Flux.  The 
curves  show  the  variation  in  response  to 
a  beam  1.6  mm  (l/l6  inch)  in  diameter  of 
unchopped  incident  flux  from  a  constant 
tungsten  source,  as  the  beam  is  scanned 
across  the  detector.  The  lower  two  curves 
show  results  obtained  when  diffusing 
screen  No.  2  was  used. 
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to  form  the  thermopile.  The  scans  in  the  A-A'  direction  were  made  across  the  center  of  the 
sensitive  area  along  the  line  between  the  two  columns.   The  scans  in  the  B-B'  direction  were 
made  across  the  center  of  the  sensitive  area,  along  the  long  axis  of  the  two  plates  in  the 
third  row.  In  both  cases  scans  were  made  (l)  using  incident  radiation  chopped  at  13  cycles 
per  second  and  an  amplifier  tuned  to  13  c.p.s.  and  (2)  unchopped  incident  radiation  and  a 
DC  amplifier.  Results  are  shown  in  Figure  26  for  the  AC  scans,  and  in  Figure  27  for  the  DC 
scans. 

In  Figure  26  for  the  A-A'  scans  it  can  be  seen  that,  with  the  screen  holder  only,  there 
are  three  distinct  peaks,  with  some  indication  of  two  others,  corresponding  to  the  position 
of  the  five  rows  of  plates.  Diffusing  screen  no.  2  was  quite  effective,  and  gave  relatively 
uniform  response  across  the  sensitive  area.  In  the  B-B'  scan  without  the  screen  holder,  it 
can  be  seen  that  the  plate  on  the  right  had  greater  sensitivity  than  that  on  the  left.  The 
screen  holder  alone  had  very  little  effect,  as  would  be  expected.  The  no.  2  diffusing  screen 
has  practically  eliminated  the  short-range  peaks  seen  in  the  other  two  scans,  but  the  greater 
sensitivity  of  the  plate  on  the  right  is  still  evident. 

In  Figure  27,  the  DC  scans  in  the  B-B'  direction  are  somewhat  similar  to  the  equivalent 
AC  scans  for  the  screen  holder  only  and  diffusing  screen  no.  2,  but  with  no  screen  holder  the  j 
curve  is  relatively  smooth,  and  shows  that  the  plate  on  the  left  has  greater  sensitivity  than 
that  on  the  right.  The  scans  in  the  A-A'  direction  show  five  distinct  peaks  with  no  screen 
holder.  The  curve  with  the  screen  holder  only  is  qualitatively  similar  to  the  equivalent 
AC  scan,  and  relatively  uniform  response  is  obtained  with  the  no.  2  diffusing  screen. 

The  variations  in  sensitivity  over  the  sensitive  area  of  the  detectors  could  readily 
account  for  the  unsatisfactory  results  reported  previously  for  the  ellipsoidal  mirror  reflec- 
tometer.  The  incident  beam  is  concentrated  in  an  area  on  the  specimen  approximately  ?mm 
square,  while  the  radiant  energy  reflected  from  the  specimen  is  focused  into  an  image  approx- 
imately 1  cm  square. 

3.3a-3  Angular  Sensitivity  of  Detectors 

Since  the  detector  needed  to  be  equally  sensitive  to  radiation  striking  at  any  angle 
from  0  to  24°,  the  angular  sensitivity  of  the  thermopile  detector  was  measured.  To  make  this 
measurement,  the  detector  was  mounted  on  a  milling  head  with  its  sensitive  area  in  a  vertical 
plane  and  with  the  two  columns  of  plates  vertical,  in  a  position  such  that  the  center  line 
of  the  sensitive  area  coincided  with  the  vertical  axis  of  rotation  of  the  milling  head.  An 
image  of  the  exit  slit  of  the  monocnromator,  3  x  3  mm  in  size,  was  focused  on  the  center  of 
the  sensitive  area  from  a  direction  normal  to  it,  by  means  of  a  6-in.  diameter  spherical 
mirror  having  a  49-in.  radius  of  curvature.  The  cone  of  rays  thus  had  a  half -angle  width  of 
about  3i°.  The  monocnromator  was  adjusted  to  give  a  band  of  flux  centered  at  2.2  microns. 
The  response  of  the  detector  was  recorded  as  Rn  when  the  axial  ray  of  the  incident  beam  was 
normal  to  the  sensitive  area.  The  milling  head  was  then  rotated  to  give  incident  angles  of 
5,  10,  15,  20,  25,  30,  40,  50,  60,  and  70°  to  the  normal,  and  the  response  of  the  detector 
was  recorded  at  each  setting  as  fig,  9  being  the  angle  of  incidence.  The  data  were  normalized 
by  dividing  each  reading  by  the  reading  at  normal  incidence,  and  plotted  as  a  function  of 
angle  of  incidence  to  produce  the  curve  shown  in  Figure  28.  Similar  measurements  were  made 
with  a  cover  plate  0.15  in.  in  thickness  with  a  1  x  1  cm  hole,  centered  over  the  sensitive 
area.  The  entire  procedure  was  then  repeated  with  flux  of  8  microns  wavelength. 

The  experimental  curves  are  compared  in  Figure  28  with  two  computed  theoretical  curve 
The  top  curve,  in  which  Rn/Rg  =  1  at  all  angles  would  be  obtained  if  the  detector  was  equally 
sensitive  to  flux  striking  it  at  all  angles,  and  if  all  of  the  incident  flux  struck  the  sen- 
sitive area.  The  lower  curve,  in  which  Rn/Rg  —  cos  9,  would  be  obtained  if  the  detector  was 
equally  sensitive  to  flux  striking  it  at  all  angles,  and  was  completely  filled  at  normal 
incidence. 

From  the  experimental  curves  in  Figure  28  it  can  be  seen  that  the  sensitivity  increases 
slightly  from  normal  to  20°,  then  decreases.  The  increase  is  undoubtedly  due  to  the  fact 
that  as  the  illuminated  area  of  the  detector  increases,  the  more  sensitive  areas,  as  shown  in 
the  previous  section,  become  illuminated.  The  sharp  drop  in  response  beginning  at  about  30° 
is  due  to  some  of  the  flux  being  lost,  either  not  admitted  through  the  window  or  not  striking 
the  sensitive  area  if  admitted.  The  presence  of  the  cover  plate  increases  the  rate  of  fall 
off  in  this  range,  as  might  be  expected. 
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0  (ANCLE    FROM    NORMAL     OF    OETECTOR) 

Figure  28.  Angular  Sensitivity  of  the 
Thermopile  Detector  in  the  A-A'  Plane  at 
2.2  and  8.0  Microns.  The  curves  show 
the  variation  in  response  to  a  beam  of 
constant  flux,  3  mm  square,  as  the  angle 
of  incidence  is  varied  in  the  A-A1  plane. 


Q   (ANGLE    FROM    NORMAL     OF    OETECTOR) 

Figure  29.  Angular  Sensitivity  of  the 
Thermopile  Detector  in  the  B-B1  Plane  at 
2.2  and  8.0  Microns.  The  curves  show  the 
variation  in  response  to  a  beam  of  con- 
stant flux,  3  mm  square,  as  the  angle  of 
incidence  is  varied  in  the  B-B'  plane, 
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Figure  30.  Diffusing  Elbow  Designed  to  Reduce 
Areal  Sensitivity  of  Thermopile  Detector. 
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Figure  31-  Angular  Sensitivity  of 
Thermopile  Detector  with  Diffusing 
Elbow,  in  the  Plane  Perpendicular  to 
the  Diffusing  Surfaces. 
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Similar  tests  were  made  with  the  detector  mounted  with  the  five  rows  of  plates  vertical. 
The  results  shown  in  Figure  29  are  similar  to  those  in  Figure  26  except  that  the  increase  in 
signal  from  0  to  20°  was  not  observed.  This  is  undoubtedly  due  to  the  fact  that  as  the  angle 
was  increased,  more  plates  were  illuminated,  but  in  the  same ' relative  areas,  hence  the  signal 
remained  constant. 

The  curves  plotted  in  Figures  28  and  29  show  no  significant  differences  for  flux  of  2.2 
and  8  microns,  respectively.  The  fact  that  the  curves  are  nearly  flat  from  C°  to  25°  indicate 
that  variation  in  angular  sensitivity  will  be  no  problem  with  the  ellipsoidal  mirror  reflec- 
tometer  since  the  marginal  rays  strike  the  thermopile  detector  at  an  angle  of  approximately 
24°  from  the  normal. 

3.3a-4-  Flux  Averaging  Devices 

The  data  on  the  variation  in  spatial  sensitivity  of  the  detectors  referred  to  in  Sec- 
tion 3 • 3a-2 " indicate  that  a  flux-averaging  device  will  be  required  for  use  with  any  of  the 
available  large-area  detectors.  The  function  of  such  a  device  is  to  spread  the  available 
flux,  regardless  of  image  size,  shape  or  intensity  distribution,  uniformly  over  the  sensitive 
area  of  the  detector. 

The  literature  provides  at  least  two  references  to  this  problem.  The  first  is  the  work 
of  Bennett  and  Koehler  [  5  1,   who  used  a  small  integrating  sphere  to  average  out  the  signal 
over  a  photomultiplier  detector.  The  second  is  the  work  of  RohzhinC6],  who  tried  light 
ducts  and  integrating  spheres  to  average  the  signal  over  the  sensitive  area  of  a  photomulti- 
plier.  However,  these  references  offer  solutions  only  in  the  ultraviolet  and  visible  pro- 
tions  of  the  spectrum,  where  good  integrating  sphere  coatings  are  available.  In  the  infra- 
red, the  spectral  region  for  which  the  reflectometer  was  specifically  designed,  no  one  has 
yet  shown  that  satisfactory  integrating  sphere  coatings  exist  for  use  beyond  about  4  microns. 
However,  Birkebak  [7]  has  shown  that  sulfur  is  both  a  good  diffuser  and  reflector  at  2  mi- 
crons and  at  4  microns,  and  he  assumes  that  it  is  usable  out  to  10  microns.  However,  he  does 
not  mention  the  specific  form  of  sulfur  that  was  used  for  his  measurements,  or  his  method  of 
applying  it  to  the  sphere  wall.  Kronstein  et  al  [8]  report  that  mu  sulfur  is  a  good  reflec- 
tor out  to  15  microns,  and  give  a  spectral  reflectance  curve.   They  did  not,  however,  use  it 
as  a  sphere  coating. 

Any  averaging  device  reduces  the  efficiency  of  a  detector  system,  because  some  of  the 
incident  flux  is  absorbed  in  the  diffuser,  and  some  is  scattered  away  from  the  sensitive 
area  of  the  detector  and  is  lost.  Since  the  flux  avialable  for  measurement  at  wavelengths 
near  1  and  15  microns  is  already  very  low,  it  is  important  that  losses  in  the  diffuser  be 
kept  as  low  as  possible.   In  general,  the  losses  in  the  diffuser  tend  to  increase  with  an 
increase  in  its  efficiency  as  a  diffuser. 

Three  different  types  of  diffusing  devices  were  investigated.  They  are  listed  in  esti- 
mated order  of  increasing  efficiency  of  diffusion  as:   (l)  a  diffusing  screen  placed  direct- 
ly over  the  detector,  (2)  a  light  duct  with  diffusing  walls  or  a  diffusing  surface  in  the 
system,  and  (3)  an  integrating  sphere  coated  with  a  material  having  high  reflectance  in  the 
infrared  and  sufficient  diffusing  power  to  permit  it  to  be  used  as  an  averaging  device. 

The  first  device  tested  was  a  roughened  sodium  chloride  window.  The  data  in  Figures  26 
and  27  indicate  that  it  was  only  partially  effective.  Another  diffusing  screen,  which  has 
not  yet  been  tested,  can  be  made  by  sandwiching  a  thin  transmitting  layer  of  f^nall  particles 
of  any  convenient  transparent  material  between  two  polished  windows  of  sodium  chloride. 
Sodium  chloride,  mu  sulfur,  and  cesium  bromide  would  be  suitable  for  such  use. 

The  second  effort  was  to  construct  a  diffusing  elbow,  shown  in  Figure  30.  This  elbow 
essentially  eliminated  the  spatial  sensitivity  when  the  flux  was  incident  on  one  end  of  the 
elbow,  and  the  detector  was  placed  at  the  other  end.  However,  Figure  31  indicates  that  the 
sensitivity  of  the  elbow-detector  combination  varied  with  angle  of  incidence  across  the  dif- 
fusing surfaces,  while  Figure  32  shows  no  such  angular  sensitivity  in  the  plane  of  incidence 
across  the  mirror  surface.  Figure  33  is  yet  another  design,  in  which  the  walls  are  all  mir- 
rors, except  for  the  diffuse  45°  surface  used  to  reflect  the  incoming  radiation  toward  the 
detector.  One  diffusely  reflecting  surface  for  use  with  this  device  consists  of  a  series  of 
spherical  depressions  in  aluminum,  each  of  1/16  in.  radius,  spaced  0.088  in.  apart  in  a  hex- 
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agonal  close  pack  array.  Goniophotometric  reflectance  curves  for  this  surface  for  white  light 
incident  at  45°  are  shown  in  Figure  34-  This  surfc.ee  has  since  been  liquid  honed,  and  then 
gold  plated.  The  liquid  honing  gives  a  diffusing  surface  of  small  roughness,  which,  in. com- 
bination with  the  large  roughness  of  the. spherical  depressions,  should  reduce  the  height  of 
the  specular  peak  at  45°  seen  in  Figure  34.  Other  diffusing  surfaces  to  be  tested  in  this 
elbow  include  mu  sulfur,  and  a  gold-plated  wire  mesh  manufactured  by  the  General  Electric  Co". 

The  final  approach  was  that  of  developing  a  small  conventional  integrating  sphere  that 
would  be  suitable  for  use  throughout  the  1  to  15  micron  wavelength  range. 

As  preliminary  trial,  a  small  integrating  sphere  diffuser,  4-5  cm  in  diameter,  coated  on 
the  inside  with  3M  white  Velvet  coating*,  was  designed  and  built.  The  entrance  aperture  and 
viewing  port  were  located  90°  apart  on  a  great  circle.  Since  the  detector  plus  diffuser  would 
not  fit  into  the  available  space  at  the  first  focus  of  the  ellipsoid,  it  was  not  possible  to 
make  direct  measurements  of  the  incident  flux.  However,  comparison  reflectance  measurements 
of  several  specimens  having  different  geometric  distributions  of  reflected  flux  gave  results 
that  were  in  reasonable  agreement  with  those  obtained  with  other  instruments.  The  3M  coating 
is  usable  only  at  wavelengths  below  about  2u-,  but  these  preliminary  measurements  indicated 
that  the  basic  design  of  the  instrument  was  sound,  and  that  it  should  be  capable  of  giving 
accurate  results  if  a  suitable  infrared  diffuser  could  be  found. 

Three  tests  were  devised  to  evaluate  the  effectiveness  of  sphere  coatings.  The  general 
optical  system  for  these  tests  is  shown  schematically  in  Figure  35.  A  152mm  (6  in.)  diameter 
spherical  mirror  of  1.24m  (49  in.)  radius  is  used  to  focus  the  beam  from  the  exit  slit  of  the 
monochromator  onto  the  entrance  port  of  the  sphere.  The  sphere  is  mounted  in  a  milling  head, 
so  that  it  can  be  moved  20  cm  (8  in.)  in  the  x  and  y  directions,  and  rotated  36O0  in  the  x-y 
plane. 

In  the  first  test,  the  sphere  was  mounted  with  its  entrance  port  at  the  center  of  rota- 
tion of  the  milling  head,  and  the  incident  beam  was  centered  on  the  entrance  port.  The 
sphere  was  then  rotated,  and  the  response  of  the  detector  was  recorded  as  a  function  of  the 
angular  position  of  the  sphere,  measured  as  the  angle  between  the  axial  ray  of  the  incident 
beam  and  the  normal  to  the* sphere  entrance  port.  In  the  second  test,  the  sphere  was  mounted 
as  before,  but  in  this  case  it  was  moved  along  the  axial  ray  of  the  incident  beam,  and  the 
detector  response  was  recorded  as  a  function  of  the  sphere  position.  In  the  third  test  a 
modified  sphere  was  used,  in  which  the  entire  illuminated  area  of  the  sphere  wall  was  viewed 
in  one  configuration. 

Data  from  the  first  test  were  normalized  by  dividing  the  signal  from  the  detector  at 
each  angle  of  view  by  the  signal  obtained  with  the  axial  ray  normal  to  the  entrance  port  (0° 
angle  of  view).  Typical  data  are  shown  in  Figures  36  through  40. 

If  a  perfect  integrating  sphere  were  tested  in  this  way,  and  the  detector  viewed  only  a 
portion  of  the  sphere  wall,  as  illustrated  in  Figure  41,  then  the  signal  from  the  detector 
would  change  as  the  directly  illuminated  spot  moves  around  the  sphere,  by  an  amount  propor- 
tional to  the  difference  in  brightness  between  the  areas  of  the  sphere  wall  that  are  directly 
illuminated  and  those  not  directly  illuminated".  If  the  illuminated  area  is  not  viewed  by  the 
detector,  none  of  the  flux  will  be  received  by  the  detector  on  the  first  reflection.  Thus, 
for  a  surface  that  approaches  an  ideal  integrating  sphere  coating,  the  curve  of  response  as 
a  function  of  angle  should  show  two  ranges  of  constant  response  with  a  smooth  monatonic  tran- 
sition. One  range  would  represent  those  angles  at  which  the  detector  wiews  none  of  the  illu- 
minated area,  and  the  other  would  represent  those  angles  at  which  it  views  the  entire  illu- 
minated area,  and  the  transition  would  represent  those  angles  at  which  the  detector  views  an 
increasing  fraction  of  the  illuminated  area. 

Jacques  and  Kuppenheim  f 1  ]  show  that  for  high  efficiency,  a  high  reflectance,  diffusely 
reflecting  sphere  wall  and  a  small  sphere  diameter  are  required.  High  sphere  efficiency  is 
required  in  this  application,  because  the  amount  of  flux  available  for  measurement  is  near 


*  Available  from  3M  Reflective  Products  Div.,  Minn.  Mining  &  Mfg.  Co.,  2501  Hudson  Rd., 
St.  Paul  19,  Minn.  Goniophotometric  measurements  made  at  NBS  show  that  this  coating 
clo.sely  approaches  a  cosine  diffuser. 
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Figure  36.  Angular  Sensitivity  of  Thermopile 
Detector  with  4-Inch  Diameter  MgO  Coated 
Sphere,  at  1.5  and  3.0  Microns. 
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Figure  38.  Angular  Sensitivity  of 
Thermopile  Detector  with  2-Inch 
Diameter  L-Shot  Gold-Plated  Sphere 
at  2.2  Microns. 
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Figure  37.  Angular  Sensitivity  of  Thermopile 
Detector  with  4-Inch  Diameter  L-Shot  Gold- 
Plated  Sphere  at  2.2  Microns. 
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the  lower  limit  of  the  useful  range  of  the  detector,  particularly  at  the  longer  and  shorter 
wavelengths,  near  1  and  15  microns.  Polished  metals  have  high  reflectance  in  the  infrared, 
but  are  not  good  diffusers.  Roughening  the  surface  of  a  metal  decreases  reflectance,  but 
increases  diffuseness.  It  may  be  possible  to  produce  a  usable  sphere  coating  by  suitable 
contouring  of  a  metal  surface,  or  by  applying  a  vacuum-deposited  metal  over  a  contoured  sub- 
strate. Birkebak  [  7  ]  showed  sulfur  to  be  diffuse  at  2  and  4  microns,  and  Agnew  and  McQuistan 
[9]  showed  it  to  be  diffuse  to  the  flux  from  a  Globar  at  wavelengths  above  and  below  4 
microns. 

A  number  of  different  spheres  were  built  and  coated.  These  included  50mm  (2-in.)  diam- 
eter spheres  with  gold  plated  L-glass  shot*  (l)  and  S-430  shot*  (2)  surfaces,  and  hand- 
pressed  mu  sulfur  (3)  coatings,  and  100mm  (4-ih. )  diameter  spheres  with  gold  plated  L-shot 
surface  (4)  and  MgO  coating  (5).  These  spheres  were  evaluated  at  several  different  wave- 
lengths, and  with  different  detector  configurations. 

The  data  for  the  100mm  (4-in. )  diameter  MgO  coated  sphere  are  shown  in  Figure  36.  At 
1.5(1,  where  MgO  is  a  good  sphere  coating,  the  signal  was  constant  to  within  1$  with  angle 
out  to  an  angle  of  60°,  and  then  it  decreased  due  to  image  spill-over  from  the  sphere  en- 
trance port.  At  3-0  M>  ,  where  the  reflectance  of  MgO  is  low  and  there  is  a  significant  spec- 
ular component,  the  signal  varied  appreciable  with  angle,   (it  should  be  noted  that  with  the 
detector  configuration  used  in  these  tests,  the  detector  viewed  at  least  a  small  portion  of 
the  illuminated  area  of  the  sphere  at  all  angles.)  Thus,  the  curves  in  Figure  36  show  a 
marked  difference  in  the  behavior  of  MgO  at  a  wavelength  where  it  is  known  to  be  a  good 
sphere  coating,  and  at  a  second  wavelength  where  it  is  known  to  be  a  poor  integrating  sphere 
coating . 

Figures  37,  38,  and  39  show  data  for  the  gold-plated  spheres.  Figure  37  is  for  the 
100mm  (4-in.)  diameter  L-shot  sphere  at  2.2  p.,  and  indicates  that  this  combination  is  not 
useful  as  an  integrating  sphere  at  this  wavelength.  In  this  case  tn£  detector  was  in  the 
plane  of  incidence.  Figure  38  is  for  a  50mm  (2-in. )  diameter  sphere  with  a  similar  coating. 
In  this  case,  the  detector  is  centered  on  a  radius  of  the  sphere  that  was  inclined  about  20° 
above  the  plane  of  incidence,  so  that  the  illuminated  area  was  directly  below  the  detector 
at  an  incident  angle  of  about  -15 °.  From  10°  to  50°  the  respones  varies  markedly  with  angle 
of  incidence,  indicating  that  the  sphere  wall  is  specular,  and  hence  not  suitable  for  an 
integrating  sphere.  From  -25°  to  0°,  there  is  little  variation  in  respones  with  angle  of 
incidence.  At  these  angles  of  incidence  the  illumiated  area  is  near  and  just  below  the 
detector  port,  where  it  is  not  viewed  by  the  detector.  Apparently,  if  the  incident  flux  can 
be  confined  to  this  area  of  the  sphere,  it  might  be  suitable  for  use  as  a  diffuser  at  2.2u.  , 
if  not  as  an  integrating  sphere.  It  remains  to  be  determined  whether  the  curve  will  be  simi- 
lar at  longer  wavelengths,  since  the  roughness  of  the  surface  is  oh  the  order  of  30n-in. 
rms,  and  the  geometric  distribution  of  reflected  energy  might  change~wlth  wavelength. 

Figure  39  shows  results  obtained  with  the  gold-plated  S-460  shot  50mm  (2-in.)  diameter 
sphere,  at  wavelengths  of  2.2,  5,  and  8  microns.   In  this  case  the  detector  port  is  in  the 
plane  of  incidence,  and  is  diametrically  opposite  the  entrance  port.  The  curves  indicate 
that  there  is  a  large  specular  component  to  the  reflected  flux,  at  45°,  and  that  the  reflec- 
tance characteristics  do  not  change  appreciably  with  wavelength.  The  S-460  shot  surface  has 
a  roughness  of  about  150g,-in.  rms,  hence  no  effect  of  wavelength  would  be  expected  in  this 
range. 

From  the  tests  reported  above,  it  is  concluded  that  the  gold  plating  on  a  roughened  sub- 
strate does  not  appear  promising  as  a  true  integrating  sphere  coating,  but  it  appears  likely 
that  suitable  diffusers  can  be  designed  for  use  with  such  surfaces. 


*  L-shot  spheres  are  available  from  the  Pangborn  Co.   They  are  glass  spheres,  -200  +  325 
mesh  SAE.  A  surface  blasted  with  these  spheres  has  a  roughness  of  about  30|i-in.  rms.  The 
S-460  shot  is  also  available  from  the  Pangborn  Co.   These  are  steel  spheres,  10-mesh  SAE,  and 
produce  a  surface  of  about  150  pi -in.  rms  roughness. 
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Figure  40  shows  data  obtained  on  the  50ram  (2-in.)  diameter  sphere  coated  with  hand-peened 
mi  sulfur.*  These  data  show  that  the  sulfur  coating  is  highly  diffuse  at  1.5,  2.2,  5-0,  and 
10. Op.  ,  since  each  of  the  curves  has  two  flat  areas,  with  a  smooth  monotonic  transition  between 
them.  This  is  what  was  predicted  earlier  for  a  diffusely  reflecting  coating.  The  tests  indi- 
cate that  a  sulfur  coated  sphere  would  be  a  satisfactory  dif fuser  if  the  sphere  were  used 
in  a  configuration  such  that  the  illuminated  area  was  either  always  totally  viewed  or  never 
viewed  by  the  detector. 

A  second  series  of  tests  was  performed  in  an  effort  to  establish  more  directly  the  suit- 
ability of  the  various  coated  spheres  for  use  with  the  ellipsoidal  mirror  reflectometer.  In 
these  tests  the  detector  port  and  entrance  port  were  centered  on  radii  of  the  sphere  135° 
apart,  as  indicated  in  Figure  35.  The  sphere  was  mounted  in  the  milling  heat  (Figure   35) 
with  the  axial  ray  of  the  incident  beam  centered  on  and  normal  to  the  entrance  port.  The 
sphere  was  moved  along  the  axial  ray  of  the  incident  beam,  and  the  detector  response  was  re- 
corded as  a  function  of  sphere  position.  This  test  was  designed  to  evaluate  the  variation  in 
detection  response  with  size  of  the  illuminated  area  on  the  sphere  wall,  when  the  total  in- 
cident flux  is  held  constant.  Since  the  incident  beam  is  diverging  from  a  focus,  the  size  of 
the  illuminated  area  could  be  varied  from  the  mini  mum  size  when  the  beam  was  focused  on  the 
sphere  wall  [or  inside  the  sphere  for  the  100mm  (4-in. )  diameter  sphere]  to  a  ma-x-iTmiTn  size 
when  the  marginal  rays  fell  just  inside  the  entrance  port.  The  maximum  size  ratio  was  6.15  to 
1  for  the  50mm  (2-in.)  diameter  sphere  and  about  2.36  to  1  for  the  100mm  (4-in.)  diameter 
sphere.  The  measured  sphere  position  was  converted  to  the  area  of  sphere  wall  illuminated  by 
the  incident  beam  from  the  beam  geometry,  and  each  area  was  divided  by  the  cross-sectional 
area  of  the  beam  at  the  focal  point  to  obtain  the  area  ratio  for  each  position.  The  detector 
response  Q  at  each  position  was  divided  by  the  response  Qq  at  the  position  where  the  beam  was 
focused  on  the  entrance  port,  to  obtain  the  response  ratio  Q/Qq.  Response  ratio  was  then  plot- 
ted as  a  function  of  area  ratio  to  obtain  the  curves  shown  in  figure  41* 

The  foregoing  test  was  designed  to  simulate  the  conditions  that  exist  when  the  detector 
is  moved  from  the  first  to  the  second  focal  point  of  the  ellipsoidal  mirror.  At  the  first 
focal  point,  the  illuminated  area  on  the  sphere  wall  will  be  on  the  order  of  10  x  12mm,  and 
at  the  second  focal  point  it  will  be  on  the  order  of  18  x  20  mm,  hence  the  test  conditions 
fully  cover  the  range  of  areas  that  will  be  encountered  in  service. 

The  data  in  Figure  41  show  that  with  the  100  mm  (4-in.)  diameter  MgO  sphere  at  1.5U-  , 
there  was  a  spread  in  response  ratio  of  about  0.8  percent.  These  results  are  consistent  with 
data  in  Figure  13,  which  show  a  spread  in  response  of  about  one  percent  with  change  in  angle 
of  incidence.  The  spread  in  response  in  Figure  41  is  in  part  due  to  the  fact  that  the  detector 
viewed  part  of  the  illuminated  area  at  some  settings. 

The  50  mm  (2-in.)  diameter  L-shot  sphere  at  2.2  p,  showed  a  spread  in  response  of  slightly 
more  than  2  percent,  which  is  above  the  acceptable  limit.  However,  it  is  felt  that  with  care- 
ful placement  of  the  illuminated  area  and  internal  baffling  of  the  sphere,  this  spread  can  be 
reduced  to  below  one  percent,  so  that  the  gold-plated  L-shot  surface  can  be  used  in  a  diffuser, 
particularly  at  wavelengths  where  the  reflectance  of  sulfur  is  low. 

The  50  mm  (2-in.)  diameter  sulfur-coated  sphere  showed  only  a  0.5  percent  spread  in  re- 
sponse at  2.2  and  8. Op,  .  Thus,  the  sulfur-coated  sphere  should  be  satisfactory  as  a  diffuser 
jout  to  at  least  8  p  .  Because  of  the  lower  reflectance  of  sulfur  at  longer  wavelengths,  its 
iuse  beyond  10.5  p.  is  questionable..  The  net  results  of  these  tests  is  to  show  that  the  sulfur- 
[coated  sphere  is  an  acceptable  averaging  device  in  the  2  -  8p,  range,  and  that  the  L-shot  gold- 
plated  sphere  may,  with  careful  design,  be  used  in  the  10  -  15 p.  range. 

The  sulfur-coated  spheres  that  were  tested  had  been  coated  by  pressing  the  sulfur  in  place 
jby  hand.  A  new  technique  has  since  been  developed  for  application  of  the  sulfur  by  spraying 
as  an  alcohol  suspension.  It  is  believed  that  the  spheres  coated  by  the  new  techniques  have 
a  more  nearly  uniform  coating  than  those  prepared  by  hand  pressing,  and,  hence,  may  show  even 
smaller  spreads  in  response  with  variation  in  size  of  the  illuminated  area. 

I 


*  Crystex  90  brand  sulfur,  purchased  from  the  Stauffer  Chemical  Co.  This  is  90£  insoluble 
(mu)  sulfur,  99-5$  elemental  sulfur,  with  3-4  p.  particle  size.  Data  on  the  reflectance  of 
sulfur  are  given  in  reference  8. 
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The  tests  described  above  indicate  qualitatively  that  sulfur  may  also  be  useful  as  an 
integrating  sphere  coating.  With  this  "Ln  mind,  a  third  test  to  quantitatively  evaluate 
integrating  sphere  coatings  was  devised.  The  results  of  the  first  test,  shown  in  Figure  40 
for  sulfur,  in  which  the  illuminated  area  was  moved  around  a  great  circle  of  the  sphere, 
gave  a  qualitative  indication  of  the  utility  of  a  surface  in  an  integrating  sphere.  However, 
the  values  for  the  difference  in  readings  when  the  illuminated  area  was  not  viewed  directly 
by  the  detector  (area  c— d  in  Figure  11,  low  flat  portion  of  the  curves  in  Figure  40 ),  and 
when  it  was  viewed  directly  by  the  detector  (area  a-b  in  Figure  11,  high  flat  portion  of 
the  curves  in  Figure  40),  were  lower  by  a  factor  of  as  much  as  5  than  would  be  predicted 
by  theory.  Apparently,  the  experimental  set-up  did  not  fit  the  theoretical  model.  The 
flux  from  the  illuminated  area  not  directly  viewed  by  the  detector  (area  c-d  in  Figure  11) 
could  reach  the  detector  by  paths  other  than  by  being  multiply  reflected  from  the  d-a-b-c 
area  viewed  by  the  detector,  by  (a)  hitting  the  lip  of  the  detector  port  and  being  reflected 
to  the  detector,  and  (b)  being  diffused  to  the  detector  by  scratches  on  the  CsBr  window. 
The  net  effect  would  be  to  increase  the  height  of  the  low  flat  portion  of  the  curves  in 
Figure  40.  The  illumination  in  the  a-b  area  of  Figure  35  is  incident  at  large  angles  of 
incidence,  where  even  the  best  diffusers  tend  to  become  somewhat  specular.  Thus,  the  flux 
tends  to  be  specularly  reflected  around  the  sphere  wall  into  the  area  b-a  to  which  the 
detector  is  less  sensitive,  instead  of  being  diffusely  reflected  to  the  detector.  The  net 
effect  would  be  to  reduce  the  height  of  the  high  flat  portion  of  the  curves  in  Figure  40. 
These  two  effects  undoubtedly  are  responsible  for  the  low  ratio  of  the  two  signals,  as  com- 
pared to  the  theoretical  model. 

Detector  and  beam  configurations  for  a  new  test  are  given  in  Figure  42.  In  this  test  both 
beams  are  incident  from  directions  near  normal,  and  flux  from  areas  not  viewed  by  the  detect- 
or can  no  longer  reach  it  by  indirect  paths.  Spheres  of  this  configuration  are  now  being 
prepared  for  test. 

3.3b.  Theory  for  Tests  of  Integrating  Sphere  Coatings. 

Figure  43  Illustrates  the  theoretical  model  to  be  studied.  The  flux  from  the  first 
reflection  will  be  considered  separately  from  that  due  to  multiple  reflections  in  the  sphere. 
To  get  a  quantitative  value  for  the  ratio  of  the  two  flat  regions  discussed  above,  it  is 
necessary  to  account  for  the  flux  that  gets  to  the  detector  in  each  case. 

Case  1.  The  detector  does  not  view  the  first  reflection  of  the  incident  flux  (Iq). 
Then  the  flux  reflected  by  the  sphere  wall  on  the  first  reflection  is 


RS  =  Vs 


(11) 
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and  none  of  this  reflected  flux  impinges  directly  on  the  detector.  This  flux  is  then  reflect- 
ed by  the  entire  sphere  (except  for  the  portion  lost  out  the  entrance  and  exit  ports).  The 
flux  reflected  twice  by  the  sphere  wall 


1  "  o  S  A.  (12) 


j  Where  A„  is  the  total  area  of  the  sphere  (As  =  A^tR2)  and  Ash  =  As  -  Ajq_  -  Aj^  (AhI  and  Aj^ 
are  the  areas  of  the  exit  and  entrance  ports,  respectively).  Of  the  flux  R^,  only  R2  reaches 
the  detector  sensing  element, 

r         1 

!  H  - '. '.  17 ^ j£!  to) 

I  where  ADy  is  the  area  of  the  sphere  fully  viewed  by  the  detector,  (Aj)y  -  A^JAs  is  the  por- 
|  tion  of  the  twice  reflected  flux  in  the  area  viewed  by  the  detector,  and  fi  is  the  diffuse 
j  configuration  factor  from  (A^y  -  Aj£?)  to  the  sensing  element  of  the  detector. 

Similarly,  on  the  third  reflection  the  amount  of  flux  reaching  the  detector  is, 
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and  on  the  n'th  reflection,  it  is, 
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The  quantities  of  flux  that  reach  the  detector  are  easily  summed  to 
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which  is  the  flux  that  reaches  the  detector  when  the  detector  does  not  view  the  first  re- 
flection. 

Case  2.  The  detector  views  the  first  reflection.  In  this  case,  the  only  additional 
flux  impinging  on  the  detector  is  that  due  to  the  first  reflection,  which  from  equation  (ll), 
is 

Rl  -  I0  ps  f2  (17) 

where  f  is  the  configuration  factor  from  the  area  illuminated  by  the  beam  to  the  detector 
sensing  element.  Since  this  is  the  only  additional  flux  impinging  on  the  detector,  the  total 
flux  reaching  the  detector  for  this  case  is 

I2  =  Ii  +  I0  Ps  f2  (18) 
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The  ratio  of  the  flux  reaching  the  detector  for  these  two  cases  is 


\     h. 


X2  =  h  +   Xo  PS  f2  "     h   PS  f2 


(19) 


■which  upon  substitution  from  equation  (16)  becomes 

h       =       i  |  (f2/f!)  [AS  -  psASH]As  (20) 

X2  "     PS<ASH>  ^DV"^ 

Therefore,  the  ratio  of  the  detector  readings  when  it  does  and  does  not  view  the  illuminated 
area  is  known  for  a  completely  diffuse  sphere  coating.  Thus,  an  easily  performed  test  can 
illustrate  the  quality  of  a  given  sphere  coating  with  regard  to  its  reflectance-diffuseness 
combination, to  the  uniformity  of  coating  reflectance  throughout  the  sphere,  and  to  whether 
the  cavity  used  is  close  enough  to  a  sphere  in  shape. 

Data  will  be  taken  with  the  sulfur-coated  sphere  shown  in  Figure  42  at  wavelengths  from 
1.5  to  lOp,.  The  data  should  give  a  quantitative  evaluation  of  the  utility  of  sulfur  as  an 
integrating  sphere  coating  in  the  infrared. 

3.3c  Calibration  of  Reflectometer 

The  ellipsoidal  mirror  reflectometer  is  capable  of  both  absolute  and  relative  measure- 
ments of  near  normal  hemispherical  reflectance  in  the  infrared.  Figure  44  illustrates  the 
several  configurations  which  are  used  to  get  the  raw  data.  IQ  is  the  incident  flux,  which 
strikes  the  specimen  at  7°  from  the  normal.  When  the  detector  is  placed  at  the  first  focal 
point,  the  incident  flux  is  measured;  when  the  detector  is  placed  at  the  second  focal  point 
and  the  sample  at  the  first,  the  flux  reflected  by  the  specimen  is  measured.  The  ratio  of 
these  two  readings  gives  the  reflectance  of  the  specimen,  except  for  the  flux  reflected 
from  the  specimen  that  does  not  reach  the  detector.  In  order  to  attain  a  high  degree  of 
accuracy,  the  losses  must  be  accounted  for  precisely.  A  convenient  method  of  accounting  for 
such  losses  is  to  prepare  an  energy  balance  for  the  system.  The  terms  for  such  an  energy 
balance  are  defined  and  discussed  below. 

A)  Mirror  Losses.  Flux  is  absorbed  by  the  mirror;  therefore,  the  reflectance  of  the 
mirror  coating  must  be  known.  This  reflectance  will  vary  with  the  angle  of  incidence,  and, 
hence,  with  position  on  the  mirror,  from  0°  at  the  apex  to  about  35°  at  the  edge.   There 
will  also  be  some  polarization  of  the  reflected  flux,  which  may  affect  the  reading.  There 
will  be  some  loss  from  scattering,  due  to  scratches,  dust,  and  other  imperfections  on  the 
mirror. 

B)  Hole  Losses.  Some  of  the  reflected  flux  will  escape  through  the  hole  in  the  mirror, 
which  admits  the  incident  beam.  This  loss  varies  with  the  geometric  distribution  of  the 
reflected  flux,  and  there  is  no  accurate  general  method  for  estimating  its  magnitude. 

C)  Sample  Shielding.  Flux  leaving  the  sample  normal  to  its  surface  will  be  reflected 
back  to  the  sample,  and,  hence,  be  blocked  from  reaching  the  detector.  Most  of  this  will 
be  lost,  but  a  very  small  fraction  will  be  reflected  by  the  specimen  and  reach  the  detector. 

D)  Diffuser  -  Ellipsoid  Interref lections.  Some  of  the  flux  incident  on  the  diffusing 
sphere  will  be  reflected  back  out  the  entrance  of  the  sphere  to  the  ellipsoid  and  back  to 
the  sphere.  This  error  will  be  largest  when  the  detector  is  at  the  first  focal  point  of 
the  ellipsoid. 
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GENERAL  TEST   CONFIGURATIONS 

Figure  44«  Test  Configurations  Used  to 
Obtain  Data  for  Calibration  of  Reflect- 
ometer. 


Figure  45-  Projection  of  Hole,  Sample,  and 
Specularly  Reflected  Image  in  First  Focal 
Plane,  Top,  and  Shield  B,  Bottom. 


E)  Atmospheric  Absorption.  The  path  lengths  for  the  incident  and  reflected  flux  will 
be  different,  hence,  atmospheric  absorption  will  introduce  errors  unless  a  non-absorbing 
atmosphere  is  used. 

F)  Sample  Holder  Losses.  Those  parts  of  the  sample  holder  and  its  supports  not  shaded 
by  the  sample  will  shade  the  detector,  and  cause  some  loss  of  flux. 

G-)  Edge  Losses.  If  the  sample  is  not  accurately  aligned  in  the  first  focal  plane  of 
the  ellipsoid,  some  of  the  flux  reflected  by  the  specimen  will  miss  the  ellipsoidal  mirror 
and  be  lost. 

Rigorous  corrections  for  the  above  losses  are  required  for  accurate  absolute  reflect- 
ance measurements,  since,  in  general,  the  corrections  are  different  for  measurement  of 
incident  and  reflected  flux.  In  comparative  reflectance  measurements,  the  errors  tend  to 
cancel  out,  since  they  are  approximately  the  same  for  the  specimen  and  standard,  especially 
if  they  have  approximately  the  same  geometric  distributions  of  reflected  flux. 

The  following  procedure  was  used  to  measure  absolute  reflectance.  The  configurations 
used  are  shown  as  parts  1,  2,  and  3  of  Figure  44»  In  Figure  kh-1,   Rq  is  measured.  Shield 
A  is  placed  above  the  sphere  entrance  port  to  admit  only  I0  and  to  eliminate  multiple  re- 
flections between  the  detector  and  the  ellipsoidal  mirror.  It  is  a  flat  plate,  painted 
on  both  sides  with  Parson's  black,  mounted  1  inch  above  the  sphere,  with  a  hole  just  larger 
than  the  incident  beam.  With  the  shield  in  place,  the  detector  reading  is  an  accurate 
measurement  of  R0  (assuming  that  the  diffuser  has  eliminated  the  effect  of  variations  in 


spatial  sensitivity  of  the  detector), 
percent . 


Use  of  the  shield  reduced  the  R0  value  by  about  4«S 
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For  the  reflected  beam,  much  more  information  is  needed.  First,  a  measurement  is  made 
with  only  the  sample  in  the  first  focal  plane,  as  in  Figure  44-2.  The  measured  flux  is  then 
a  rough  estimate  of  the  reflected  flux,  which  must  be  corrected  for  mirror  losses,  hole  losses, 
sample  shading,  atmospheric  absorption,  and  interreflections  between  the  diffuser  and  the 
ellipsoidal  mirror.  The  arrangement  shown  in  Figure  44-3  is  used  to  obtain  additional  informa- 
tion about  the  flux  incident  on  the  ellipsoidal  mirror  near  the  hole.  There  is  a  one-to-one 
correspondence  between  areas  on  the  ellipse  and  areas  in  the  first  focal  plane.  Figure  45 
shows  the  projection  of  the  hole,  sample,  and  specularly  reflected  image  in  the  first  focal 
plane.  The  hole  loss  and  specimen  shielding  loss  can  be  computed  from  the  flux  incident  on 
these  areas.   Shield  B  (Figure  45)  is  inserted  into  the  first  focal  plane,  and  the  detector 
reading  is  taken  as  R2.  The  reading  without  any  shield  is  Ri.  Then  R3  =  R]_  -  R2  is  a  meas- 
ure of  the  flux  striking  shield  B.  If  we  assume  that  the  flux  is  uniformly  distributed  over 
the  solid  angle  represented  by  the  shield,  or  varies  linearly  over  the  solid  angle,  we  can 
correct  for  hole  and  sample  shielding  losses  as  follows.  The  total  flux  striking  the  shield, 
of  area  Ajljg,  is  Ro,  but  the  flux  that  would  have  struck  the  area  Aj!j,  the  projected  area  of 
the  hole,  has  passed  out  the  hole.  Thus,  the  flux  density,  assuming  uniform  distribution,  is 
R3/AHS  ~  A^,  and  the  corrected  R3  would  be  (R3Ahs)/(Ahs  ~  Ah).  The  sample  also  shades  the 
detector  with  an  area  A's,  which,  if  we  assume  it  to  have  the  same  flux  density  incident  on  it 
from  the  ellipsoidal  mirror,  will  increase  the  corrected  R3  to         A'  +  A1 


3  ^S  "  M 


The  corrected  reading,  R  will  then  be 


R  =  R„  +  R„ 
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+  a;. 
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Some  of  the  flux  that  strikes  the  sample  a  second  time  will  be  reflected  to  the  ellipsoidal 
mirror  and  the  detector,  and  thus  increase  the  R]_  and  R2  readings.  However,  this  correction 
is  considered  at  this  time  to  be  negligibly  small.   It  will  be  accurately  evaluated  in  the 
near  future. 

The  signals  read  by  the  detector  in  the  configurations  shown  in  Figure  44-2  and  44-3 
will  be  high,  due  to  interreflection  between  the  diffuser,  sphere,  and  specimen.  Computations 
indicate  that  the  error  should  be  on  the  order  of  2  to  3  percent.  Since  there  does  not  appear 
to  be  any  convenient  way  to  eliminate  the  error,  it  will  be  necessary  to  evaluate  the  error 
experimentally  and  correct  for  it. 

No  absolute  measurements  are  being  made  at  this  time  in  the  wavelength  ranges  where  at- 
mospheric absorption  is  significant.  Where  the  error  is  significant,  absolute  measurements 
will  be  made  in  a  non-absorbing  atmosphere.  The  error  cancels  out  in  comparison  measurements, 
hence  a  standard  that  has  been  calibrated  in  a  non-absorbing  atmosphere  will  be  used  for  such 

measurements. 

The  remaining  major  correction  is  for  mirror  losses.  The  coating  on  the  ellipsoidal  mir- 
ror is  aluminum,  applied  in  one  second  at  a  pressure  of  10~5  torr  to  a  thickness  sufficient 
to  be  opaque.  There  is  no  overcoat.  The  reflectivity  thus  should  approach  that  of  pure  evap- 
orated aluminum,  according  to  Bennett,  et  al  [10],  and  Hass  [11].  However,  there  is  some 
scatter,  due  to  imperfections  in  the  finish  on  the  mirror,  which  will  increase  the  apparent 
absorption.  Hence  it  will  be  necessary  to  measure  the  effective  reflectance  as  a  function  of 
location  on  the  mirror.  This  will  be  done  as  follows:  a  gold  mirror  and  an  aluminum  mirror, 
of  known  reflectance  and  minimum  scatter,  will  be  used  as  samples.  They  will  be  mounted  on 
a  specimen  holder  that  can  be  tilted  to  any  angle  in  any  direction,  so  that  the  specularly 
reflected  beam  can  be  directed  to  any  part  of  the  ellipsoidal  mirror.  Measurements  of  the 
reflectance  of  the  samples  will  thus  give  the  data  required  to  evaluate  the  reflectance  of  the 
mirror  in  the  illuminated  area.  Measurements  will  be  made  at  all  wavelengths  of  interest. 
If  significant  variation  in  reflectance  is  found  over  the  surface  of  the  mirror,  it  can  be 
corrected  for.  Effects  of  polarization  by  the  ellipsoidal  mirror  will  also  be  considered. 

The  losses  due  to  shading  by  the  sample  holder  and  support  are  believed  at  this  time  to 
be  negligibly  small.  This  will  be  verified  experimentally. 
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All  of  the  errors  discussed  above  will  be  present  when  the  instrument  is  used  for  absolute 
measurements  -  i.e.,  when  the  incident  flux  is  measured  directly.  If  comparative  measurements 
are  made  -  i.e.,  the  incident  flux  is  evaluated  by  measuring  the  reflectance  of  a  standard  of 
known  reflectance,  as  indicated  in  Figure  44-2,  3>  k,   and  5  -  the  corrections  tend  to  cancel 
out  unless  there  is  a  large  difference  in  the  geometric  distribution  of  the  flux  reflected  by 
the  sample  and  standard. 


4.   Summary  of  Work  on  Relation  of  Thermal 
Emittance  to  Other  Properties 

4.1  Introduction 

One  of  the  long-range  goals  of  this  program  from  the  beginning  has  been  to  obtain  a  bet- 
ter understanding  of  the  fundamental  processes  involved  in  the -interaction  of  electromagnetic 
radiation  with  matter.  Attainment  of  such  a  goal  is  intrinsically  desirable;  further,  it  can 
contribute  to  the  attainment  of  other  more  immediate  objectives  in  two  ways.  First,  such  a 
knowledge  can  aid  in  the  development  and  evaluation  of  equipment  and  techniques  for  measuring 
the  thermal  radiation  properties  of  materials.   Second,  if  meaningful  equations  can  be  found 
or  derived  that  relate  thermal  radiation  properties  to  other  more  easily  evaluated  properties 
of  materials,  they  can  be  used  to  interpolate  and  extrapolate  thermal  radiation  property  data, 
obtained  at  a  relatively  few  temperatures  and  wavelengths,  over  much  wider  temperature  and 
wavelength  ranges. 

Pertinent  information  on  this  subject  has  been  published  by  many  authors,  representing 
widely  different  scientific  disciplines.  However,  there  does  not  appear  to  be  any  single 
compilation  of  the  major  contributions.  This  summary  includes  a  review  of  the  pertinent  lit- 
erature, including  a  discussion  of  the  limitations  of  the  various  theoretical  treatments.  In 
addition,  the  optical  behavior  of  rhodium,  for  which  reliable  data  are  available,  is  treated 
in  some  detail. 

4.2  Consideration  of  Quantum  Theory  [12] 

Quantum  theory  is  used  to  explain  on  an  atomic  or  sub-atomic  scale  the  generation  and 
absorption  of  electromagnetic  radiation.  According  to  this  theory,  electromagnetic  radiation 
consists  of  discrete  particles  called  photons,  each  of  which  has  a  given  kinetic  energy,  which 
is  related  to  frequency.  These  photons  are  generated  and  absorbed  as  units,  by  statistical 
processes.  They  travel  in  straight  lines,  at  a  velocity,  c,  which  in  vacuum  is  a  fundamental 
natural  constant,  the  speed  of  light.  Their  direction  of  propagation  is  changed  in  a  predict- 
able manner  in  passing  through  an  interface  separating  two  volumes  having  different  optical 
properties . 

Because  the  generation  and  absorption  of  photons  are  statistical  processes,  it  is  not 
possible  to  predict  when  or  where  either  process  will  occur,  but  it  is  possible  to  predict 
statistically  what  fraction  of  the  photons  passing  through  a  material  will  be  absorbed  in 
unit  path  length,  and  how  many  photons  of  a  given  energy  will  be  generated  per  unit  volume 
and  unit  time  for  any  given  material  under  given  conditions. 

Quantum  theory  has  been  used  successfully  to  account  quantitatively  for  the  emission  and 
absorption  of  gases,  and  qualitatively  to  account  for  certain  features  of  the  spectra  of  liq- 
uids and  solids.  In  general,  however,  the  mathematical  treatment  is  too  cumbersome  to  permit 
computation  of  the  complete  spectra  of  solids  and  liquids,  even  by  use  of  modern  computer  tech- 
niques . 

Classical  electromagnetic  wave  theory,  on  the  other  hand,  can  be  used  to  explain  the 
radiant  heat  transfer  on  a  macroscopic  scale.  While  it  does  not  explain  the  specific  mechan- 
isms involved  in  the  generation  and  absorption  of  individual  photons,  it  does  account  quanti- 
tavely  for  the  net  exchange  of  energy  by  radiation.  Hence  the  discussion  that  follows  will 
be  based  primarily  on  classical  electromagnetic  wave  theory. 
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4.3  Electromagnetic  Theory  of  Radiation 

Every  object  at  a  temperature  above  0  °K  is  continually  emitting  radiant  energy.  This 
energy  leaves  the  object  in  the  form  of  waves,  which  are  characterized  by  an  electric  field, 
a  magnetic  field,  polarization,  and  frequency.  The  electric  and  magnetic  fields  are  at  right 
angles  to  each  other  and  to  the  direction  of  propagation.  Polarization  refers  to  the  orienta- 
tion of  the  electric  field;  in  a  plane  polarized  beam  of  light  the  electric  fields  of  all  the 
waves  are  parallel.  The  wave  travels  in  free  space  at  the  velocity  of  light,  c.  The  nature 
of  the  wave  is  mathematically  described  by  the  Maxwell  equations,  which  can  be  found  in  any 
standard  text  on  light. 

Let  us  consider  a  plane  interface  separating  free  space  from  a  volume  consisting  of 
some  optically  homogeneous  material  as  yet  unspecified.  As  an  electromagnetic  wave  impinges 
on  this  interface  from  free  space,  interaction  of  the  wave  and  material  will  occur,  and  the 
character  of  the  wave  will  change.  It  will  be  separated  into  two  components,  one  reflected 
by  the  interface  and  one  continuing  into  the  material.  We  need  not  be  concerned  with  the 
mechanism  of  the  change,  but  only  to  quantitatively  describe  it  and  relate  it  to  the  proper- 
ties of  the  material. 

The  change  in  the  continuing  component  is  most  simply  described  in  terms  of  phase  and 
velocity  of  propagation.  The  wave  within  the  material  is  described  by 

E  =  Eoexp(j(jut  -  yx)  (22) 

in  which  E  is  the  amplitude  of  the  electric  field  at  time  t,  EQ  is  the  amplitude  at  zero  time, 
the  term  exp  (  )  indicates  that  e,  the  base  of  natural  logarithms,  is  raised  to  th3  indicated 
power,  in  parenthesis,  j  =/  -1,  w  is  frequency,  Y  =  <*  +  j3,  where  or  is  the  attenuation  factor 
and  3  is  the  phase  factor,  andx  is  the  distance  traversed  in  time  t.  If  zero  time  is  taken 
as  the  time  the  wave  impinges  on  the  interface,  E0  is  the  value  for  the  wave  in  free  space. 

The  reflected  energy  can  be  computed  by  using  boundary  conditions  compatible  with  con- 
tinuity of  the  field  vector  components.  This  approach  develops  the  Fresnel  equations  for 
the  two  reflected  components  that  are  polarized  normal  and  parallel  to  the  plane  of  incidence 
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m   pcos  Y  -Z  cos  qn 
'   n    cos  Y  +Z  cos  cpj 

«/*   v\    ."cos  cp  -Z  cos  Y  I2  ,„,  N 

P<V^>p  -  !_co.  I   +Z  cos  Y_i  (24) 


in  which  p  =  reflectance,  the  subscripts  ^ ,  V    and  \  denote  angle  of  incidence,  angle  of 
refraction  and  wavelength,  respectively,  and  the  subscripts  n  and  p  denote  the  components 
polarized  normal  and  parallel  to  the  plane  of  incidence.  Z  is  the  impedance  ratio  of  the 
material,  Z  =  Z*/Z0*,   where  Z-»-  is  the  intrinsic  impedance  of  the  material  and  Z0-k-  is  that  of 
free  space.  The  angles,  as  are  all  angles  in  optics  and  wave  theory,  are  measured  from  the 
normal  to  the  surface  to  the  direction  of  propagation.  The  direction  of  propagation  of  the 
reflected  ray  is  always  in  the  plane  of  incidence,  and  the  angle  of  reflection  is  equal  to 
the  angle  of  incidence,  but  lies  on  the  side  of  the  normal  opposite  the  incident  ray.  The 
angles  Y  and  co  are  related  by  Snell's  law 

sin  cp  =  |z|  sin  Y  (25) 

It  is  evident  that  the  reflectivities  given  by  equations  (23)  and  (24)  are  not  complex 
quantities,  but  are  real  numbers,  and  as  such  the  value  of  the  reflectivity  at  a  single 
wavelength  and  angle  of  incidence  cannot  completely  describe  the  reflection  process.  To 
evaluate  Z  or  Z#  from  measured  reflectivities  requires  information  at  various  angles.  Even 
then,  the  computation  requires  a  tedious  trial  and  error  approach  and  precise  information 
on  the  angular  dependence  of  the  reflectivity  data. 
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Further  manipulation  of  wave  relations  and  the  Maxwell  equations  involves  complex  pa- 
rameters (with  real  and  imaginary  components)  that  describe  the  wave  after  traversing  a  known 
path  length  in  the  material.  A  list  of  parameters  that  have  been  used  for  this  purpose  is 
given  in  Table  5,  with  a  brief  explanation  of  their  physical  significance. 

To  obtain  emissivity  values  for  opaque  materials,  Kirchhoff ' s  law  is  applied  to  give 

e(cp,\)  -  1  -  p(cp,\)  (26) 

where  e  (<p,  \)   is  the  spectral  emittance  at  wavelength  \  and  direction^,  and  P  (<P,\)  is  the 
spectral  specular  reflectance  for  radiation  incident  at  angel <P. 

A  special  case  of  great  interest  (and  less  complexity)  is  when  cp  =  0,  giving 

cn(0,X)  -  ep(0,\)  .  1  -  [f~f]  (27) 

Equation  (27)  can  also  be  written  in  terms  of  indices  of  refraction  and  absorption 

where  n  is  index  of  refraction,  K  is  absorption  index,  and  n*  is  the  complex  index  of  re- 
fraction, n*  =  n(l-jk).  Typical  values  of  these  parameters  are  given  in  Table  6. 

So  far,  no  assumptions  have  been  made  as  to  the  nature  of  the  material,  except  that  it 
|  is  optically  homogeneous  and  isotropic,  with  a  perfectly  smooth  interface.* 

4.4  Metals  [12,  13,  14] 

In  the  preceding  section  of  treatment  of  the  relation  of  thermal  radiation  properties  of 
materials  was  developed  on  the  basis  of  complex  factors,  without  consideration  of  the  nature 
of  the  material.  To  evaluate  any  of  the  parameters  given  in  Table  5,  it  is  necessary  to 
examine  in  detail  the  effect  of  the  impinging  electromagnetic  wave  upon  the  electron  configu- 
ration of  the  material.  The  following  discussion  is  confined  to  metallic  conductors,  and 
is  limited  in  scope  to  theories  within  the  framework  of  classical  electron  physics. 

4.4a  Free  Electron  Models 

The  earliest  attempts  to  predict  the  optical  properties  of  metals  were  made  by  Lorenz 
[15],  Drude  [16],  Kronig  [17],  and  Mott  and  Zener  [18 J,  who  assumed  the  metal  to  contain 
electrons  which  were  essentially  free  to  move  under  the  influence  of  the  electric  field  in- 
duced by  the  incident  electromagnetic  wave.  These  free  electrons  are  the  valence  electrons 
in  the  outer  shell  of  the  atoms  constituting  the  metal.  When  the  wave  is  incident  upon  its 
surface,  an  oscillating  electric  field  parallel  to  the  surface  is  induced  in  the  metal.  The 
free  electrons  will  oscillate  under  the  influence  of  this  field  at  the  frequency  of  the  in- 
cident wave.  There  is  a  phase  difference  between  the  oscillation  of  the  electrons  and  that 
! of  the  field,  caused  by  a  viscous  damping  force,  Yf  ,  arising  from  the  collisions  between 
accelerated  electrons  and  the  atomic  lattice. 

The  equation  of  motion  for  an  electron  moving  under  these  conditions  is: 

my"  +  2nmYfy'  =  q  Eexp  (  jut)  (29) 

where  m  is  the  mass  of  an  electron,  y'  is  velocity,  y"  is  the  acceleration,  q0  is  the  charge 
on  an  electron,  E  is  the  amplitude  of  the  electric  field,  exp  (__)  indicates  that  the  base  of 
natural  logarithms  is  raised  to  the  power  in  parentheses,  j  =/  -1,  w   is  frequency  and  t  is 
time.  The  current  density  corresponding  to  this  electronic  motion  can  then  be  computed, 
which  gives  the  relationship  between  the  parameters  descriptive  of  the  wave  and  the  parame- 
ters descriptive  of  the  electronic  configuration.   The  most  popular  form  is  in  terms  of  the 


-;;-  It  is  also  assumed  that  the  magnetic  effects  are  negligible  compared  to  the  electric 
effects,  and  that  the  permeability,  M-*,  of  the  material  is  the  -same  as  that  of  free  space. 
These  are  reasonable  assumptions  for  most  materials  of  primary  concern  here. 
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dielectric  constant,  K* 


1  +  J7 
K*  =  1  -  (\/\  )3  &_  „  (30) 
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where  X^  wavelength  of  incident  wave,  \Q   =  c/qg  (jjr)  ,  c  =  speed  of  light,  Nf  =  number  density 
of  free  electrons  per  unit  volume,  >f  =  2tjct,  =  c/Xf,  t  =  relaxation  time. 

In  order  to  predict  the  optical  properties,  two  parameters  of  the  material  must  be  eval- 
uated, t  and  Nf .  If  Nf ,  the  number  of  free  electrons  excited  by  the  induced  field,  is  con- 
sidered to  be  the  number  of  valence  electrons  per  unit  volume,  the  relaxation  time,  T,  can  be 
computed  from  electrical  conductivity  and  the  assumption  of  a  spherical  Fermi  surface.  This 
approach  gives  favorable  results  compared  with  measured  values  for  liquid  metals,  but  fails 
for  more  complex  metals.  Figure  46  shows  the  general  shape  of  the  computed  and  measured 
emissivity  and  optical  constants  for  mercury,  for  which  the  theory  gives  fair  agreement. 

If  it  is  assumed  that  the  phase  change  arising  from  electronic  collisions  can  be  neglect- 
ed, the  theoretical  model  is  greatly  simplified,  and  the  optical  constants  and  the  emissivity 
can  be  expressed  in  terms  of  the  d.c.  electrical  resistivity  alone.  Such  an  assumption  was 
made  in  the  Drude  model  referred  to  above  at  very  long  wavelengths,  or  for  conditions  where 
the  relaxation  time  is  infinite.  The  resulting  equation  for  these  conditions  is  the  celebrated 
Hagen-Rubens  Tl9]  (sometimes  called  the  Drude  Zener)  relation: 

«(0,\)  ,  0.365  (Pm/>0   -  0.0464  (p  /\)  (31) 

where  e(0,X)  is  the  normal  speGtral  emissivity  andpm  is  the  d.c.  electrical  resistivity.  The 
usefulness  of  equation  (31)  has  long  been  recognized,  but  its  limitations  have  not  always  been 
appreciated.  Comparison  of  values  computed  from  this  equation  with  experimental  data  reveals 
differences  which  vary  with  the  metal  involved,  the  wavelength  and  the  temperature.  These 
differences  arise  because  the  assumptions  made  in  deriving  the  equation  are  not  completely 
valid,  particularly  at  short  wavelengths. 

Many  other  equations  have  been  derived  on  the  basis  of  equation  (31),  for  computation  of 
hemispherical,  total  and  directional  emissivity.  Table  7  lists  the  equations,  together  with 
a  brief  explanation  of  their  limitations. 

4.4b  The  Multi-Free/Bound  Model 

The  failure  of  the  free  electron  model  to  yield  computed  values  that  agree  with  experi- 
mental data  in  the  wavelength  regions  of  interest  should  not  be  interpreted  as  a  failure  of 
the  classical  approach,  but  rather  should  be  ascribed  to  the  over-simplification  of  the  mech- 
anism of  the  reflection  (or  absorption)  process.  Quantum  theory  indicates  that  both  free  and 
bound  electrons  exist  in  a  metal,  and  both  types  of  electrons  exist  in  different  energy  states 
As  such,  they  react  in  different  degrees  to  the  influence  of  the  electric  field  induced  by  the 
incident  radiation.  The  free  electrons  of  a  particular  energy  state  are  characterized  by  a 
particular  viscous  damping  parameter,  Xf  (or  relaxation  time  T),  and  number  density,  Nf. 

Contributions  due  to  the  influence  of  the  induced  electric  field  on  bound  electrons  must 
also  be  considered.  They  can  be  assumed  to  be  elastically  bound  in  their  equilibrium  position; 
and  each  type  of  bound  electrons  is  characterized  by  particular  values  of  three  parameters  - 
elastic  restoring  force,  KD,  viscous  damping  force  XD,  and  number  density,  ND.  Then  a  clas- 
sical model  for  the  most  general  case,  where  contributions  due  to  several  types  of  free  and 
bound  electrons  must  be  considered,  will  have  the  following  equations  of  motion. 

Free  electrons,  type  g    g  =  1,  2,  3,  •  •  • 


my 
g 


+  2nm\gys"   qoKo"P^»t5  (32) 
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COMPARISON    OF    THE    REFLECTIVITY   OF    RHODIUM 


Figure  47  •  Comparison  of  the  Reflectivity  of  Rhodium  with 
Values  Predicted  by  the  Multi-Free/Bound  Electron  Model. 
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Figure  48.  The  Effect  of  Individual  Terms  of  the  Multi-Free/ 
Bound  Electron  Model  on  Emissivity. 


Figure  46.   (a)  Experimental 
Values  of  n  and  k  for  Mer- 
cury [45,  46,  47,  48,  49] 
Compared  with  Values  Com- 
puted by  the  Drude  [16]  Theory. 

(b)  Experimental  Values  of  Reflectivity  of  Mercury-Glass  and  Mercury-Quartz 
Interfaces  (Points)  Compared  with  That  Computed  from  the  Drude  Theory  (Curves). 

(c)  Experimental  Values  of  Reflectivity  of  Mercury  -  NaCl  Interface  (Points) 
Compared  with  That  Computed  from  the  Drude  Theory  (Curve).  Figures  from 
Reference  50. 
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h       \ . 

Here,  7  is  electron  displacement.  The  development  of  the  optical  properties  from  the  equations 
of  motion  proceeds  as  described  above  for  the  simpler  free  electron  model.  The  resulting  form 
for  the  complex  dielectric  constant,  K*,  is 


K»  =  1  + 


,2  a 


4TT8 


„   3 

me  c 


NbhN>h 

L              oh          \>, 

(34) 


where  X    =  viscous  damping  parameter  for  free  electrons,  type  g 
X , ,   =  viscous  damping  parameter  for  bound  electrons,  type  h 

K^h  =  elastic  restoring  force  on  bound  electrons,  type  h 

.  m2 

oh  =  2iic  r   ,  resonant  wavelength  of  bound  electrons,  type  h,  with  no  damping  forces  : 
^bh 

present. 
e  o   =  permittivity  of  free  space 

In  order  to  evaluate  the  optical  properties  by  means  of  these  cumbersome  equations,  a 
method  for  evaluating  Xf,  X,  ,  K^  and*  0  must  be  found.   So  far  no  theory  has  been  developed 
for  estimating  these  parameters  from  the  electronic  configuration  of  the  metal.  Such  an  ap- 
proach would  require  quantum  mechanical  treatment,  or  comparison  with  experimental  data. 

Roberts  L  20,  21,  22J ,  using  the  experimental  data  approach,  successfully  determined 
parameters  for  a  number  of  metals.  The  form  of  the  equation  is  somewhat  different  than  that 
of  equation  (34)  >  but  a  comparison  of  the  forms  and  meanings  of  equivalent  parameters  is  easily 
seen  in  Tables  8  and  9.  We  will  not  discuss  the  method  used  to  evaluate  the  parameters  from 
experimental  data,  but  the  results  are  shown  in  Table  10. 

The  success  of  the  Roberts  treatment  in  producing  computed  emissivity  values  that  agree 
in  detail  with  experimental  data  has  been  substantial.  Inflection  points  and  "humps"  have 
been  properly  described,  so  that  if  the  material  parameters  are  known,  optical  behavior  is 
completely  described  in  the  wavelength  region  considered.  No  equations  have  yet  been  developed" 
on  the  basis  of  the  multi-free/bound  electron  model  that  correspond  to  those  developed  for"  the 
simpler  free  electron  model.  It  is  apparent  that  until  equations  are  developed  that  include 
the  contributions  from  resonance  effects,  accurate  prediction  of  emissivity  from  structural 
factors  (or  good  agreement  between  computed  and  experimental  data)  is  not  to  be  expected,  par- 
ticularly at  the  shorter  wavelengths  (higher  frequencies). 

4.4c  The  Kramer-Kronig  Relation  [  23  3 

It  was  stated  above  that  two  parameters  are  required  to  completely  describe  the  optical 
behavior  of  a  material  at  a  given  wavelength.  The  Kramer-Kronig  relation  provides  an  alter- 
nate method  of  obtaining  full  knowledge  of  the  optical  behavior  from  knowledge  of  only  one 
parameter  (such  as  the  normal,  spectral  reflectivity)  at  all  wavelengths,  or,  more  precisely, 
over  the  wavelength  range  where  variations  in  the  property  occur.  This  technique  is  a  very 
powerful  tool  in  solid  state  physics  studies,  and  is  amply  described  and  demonstrated  in  the 
literature  [  24,  25 1.  A  brief  derivation  is  given  in  Appendix  A. 
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The  Kramer-Kronig  equation  has  the  form: 

1  "  v  (X) 
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where  9(X0)  is  the  phase  component  of  the  complex  quantity,  V,  which  has  the  real  part  v(X0) : 

I  T -*■»<*•>-  i£fi  (36) 

I  and  from  previous  relations  (eq.  28)  it  is  apparent  that  p(0,X)  =  V2.  (37) 

|  To  compute  the  phase  9(XQ),at  X  =  X0,  equation  (35)  is  integrated  over  the  entire  wavelength 
I  region  (of  consequence)  about  the  pole  X  =  XQ.  Knowing  the  phase  it  is  then  possible  to  solve 
for  the  real  and  imaginary  parts  of  n*. 

The  greatest  difficulty  encountered  in  using  this  technique  is  to  establish  the  wavelength 
range  over  which  the  integral  must  be  evaluated.  At  first  thought  one  might  expect  that  to 
evaluate  the  phase  at  some  longer  wavelength,  such  as  lOp.,  contributions  to  the  integral  within 
the  ultraviolet  region  may  be  ignored.  In  general  this  is  not  true,  as  can  be  verified  upon 
closer  inspection  of  the  form  of  the  integral.  As  a  rule,  full  knowledge  of  the  reflectivity 
spectrum  is  required.  The  lower  limit  of  the  significant  range  is  reached  when  the  reflectiv- 
ity goes  monotonically  to  a  zero  value  as  the  wavelength  approaches  zero.  The  upper  limit  is 
reached  when  the  reflectivity  becomes  constant  or  approaches  unity  as  the  wavelength  increases. 
The  important  consequence  of  this  is  that  the  integral  must  be  evaluated  at  «.n  wavelengths 
(especially  in  the  short-wave  region)  where  optical  activity  occurs  such  as  absorption  edges 
or  resonance  effects. 

4.5  Optical  Behavior  of  Rhodium 

Very  early  in  this  study,  data  on  the  normal  spectral  reflectivity  of  rhodium  were  col- 
lected [26].  The  choice  of  rhodium  was  based  upon  the  availability  of  reliable  data  over  a 
wide  wavelength  range.  The  reflectiveity  spectrum  from  0.058m-  to  10m.  was  felt  to  be  adequate 
information,  and  also  showed  interesting  inflection  points  in  the  wavelength  region  of  signifi- 
cance from  a  thermal  radiation  standpoint. 

Using  the  multi-free/bound  electron  model  similar  to  that  of  Roberts,  attempts  were  made, 
to  determine  the  values  of  the  parameters  which  would  give  a  reflectivity  spectrum,  figure  47, 
comparing  favorably  with  the  rhodium  data.  In  attempting  to  fit  the  observed  values  with  a 
formula  containing  just  one  free-electron  term,  the  coefficients  of  that  free  electron  term 
were  determined  first  to  produce  a  good  fit  at  the  upper  end  of  the  data —  near  X  =  10m-  — 
and  in  the  neighborhood  of  X  =  0.25m-,  where  the  reflectivity  rises  sharply.  Then  bound- 
electron  terms  were  introduced,  one  at  a  time,  to  allow  fitting  of  the  "humps"  in  the  measured- 
reflectivity  curve,  such  as  the  curve  near  0.45  m-  Each  time  a  new  term  was  introduced,  the 
parameters  of  the  older  terms  had  to  be  modified  so  as  to  maintain  the  previously  achieved 
fit;  this  in  turn  necessitated  modifying  the  parameters  of  the  new  term,  which  led  to  further 
modifications  of  the  old  terms,  etc.  It  was  originally  thought  that  the  effect  of  each  bound- 
electron  term  would  be  quite  local*,  so  that  these  modifications  would  be  small,  and  after 
one  or  two  modification  cycles  a  good  fit  would  be  obtained  both  at  the  newly-considered  "hump" 
and  at  all  locations  previously  considered.  This  was  only  partly  borne  out;  a  bound-electron 
term  introduced  to  account  for  an  apparent  slight  hump  between  li*  and  2M-had  large  effects  as 


*  It  can  be  easily  shown  that  the  effect  of  a  resonance  about  a  wavelength  X  =  X  s^  is 
not  necessarily  localized,  but  the  influence  may  extend  to  other  portions  of  the  spectrum. 
An  example  of  this  is  shown  in  Figure  48  where  the  resulting  emissivity  is  shown  including 
successively  more  terms.  Close  inspection  reveals  how  the  influence  of  an  additional  bound 
term  gives  rise  to  a  widespread  effect.  This  particular  example  is  for  tungsten  at  1800  °K 
using  the  parameters  determined  by  Roberts. 
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far  down  as  0.l6u.  The  comparison  of  various  attempts  with  the  experimental  data  is  shown  in 
Figure  47  and  numerical  values  of  the  parameters  are  recorded  in  Table  11. 

It  is  seen  from  Figure  47  that  curves  1,  2,  and  3  compare  favorably  with  the  data  except  [ 
at  wavelengths  near  and  below  0.1^.  The  optical  constants  computed  from  the  parameters  in  the, 
visible  portion  of  the  spectrum  (where  the  computed  curve  shows  good  agreement  with  experiment- 
al data)  do  not  agree  with  measured  optical  constants  for  the  same  wavelengths.  Hence  it  is 
evident  that  the  parameters  are  mathematically  expedient,  and  have  no  physical  significance. 

The  Kramer-Kronig  relations  can  be  used  to  obtain  the  phase  and  the  optical  constants  of 
a  material  at  any  particular  wavelength  from  the  normal  reflectivity  of  the  meaterial  over 
the  full  wavelength  spectrum.  When  the  observed  values  of  V  =  /  p  ( 0,X )  for  the  wavelength 
range  0.14  to  10. Op  were  used  with  equation  (36)  to  compute  the  phase  and  optical  properties, 
unrealistic  values  were  obtained.  The  computed  value  of  6(^)  at  0. 546 p,  was  0.274,  while  the 
experimentally  determined  value  is  0.395.  The  discrepancy  is  due  to  the  absence  of  data  on 
V  at  wavelengths  below  0.14M-.  Rough  calculations  using  assumed  shapes  for  the  reflectivity 
curve  below  0.14m-  showed  that  the  contribution  of  this  portion  of  the  spectrum  to  the  value 
of  9  (X)   is  apt  to  be  of  the  same  order  of  magnitude  as  the  contributions  of  the  entire  spec- 
trum at  wavelengths  beyond  0.14^. 

This  observation  further  implies  that  the  attempt  to  fit  the  reflectivity  curve  by  use 
of  the  general  classical  model  could  not  fully  succeed  because  of  the  lack  of  reflectivity 
data  at  very  short  wavelengths.  It  is  apparent  that  there  are  "humps"  in  the  reflectivity 
curve  at  wavelengths  below  0.141^;  the  bound-electron  terms  that  should  be  introduced  into  the 
equations  to  account  for  these  humps  would  modify  the  parameters  of  the  free-electron  terms 
previously  considered.  However,  since  the  locations  and  shapes  of  these  humps  are  not  predict-, 
able  from  the  reflectance  data  at  longer  wavelengths,  it  is  impossible  to  take  them  into  ac- 
count so  that  the  parameters  would  be  physically  significant. 

4.6  Summary  and  Conclusions 

This  brief  review  of  the  fundamentals  of  wave-matter  interaction  and  the  optical  behavior 
of  metals  has  focused  attention  upon  three  important  aspects  of  theoretical  predictions  of 
thermal  radiation  properties. 

(1)  A  complete  description  of  the  optical  behavior  requires  measurement  or  prediction 
of  two  parameters  at  any  wavelength  since  two  parameters  completely  characterize  the  form  of 
the  electromagnetic  wave  interacting  with  the  material. 

(2)  The  thermal  radiation  property  reflectivity  (or  emissivity)  is  a  thermodynamic  prop- 
erty resulting  from  energy  conservation  considerations  and  one  such  value  at  a  wavelength  is 
not  sufficient  information  to  completely  describe  optical  behavior.  t 

(3)  The  Kramer-Kronig  relation  provides  a  means  of  determining  optical  behavior  from 
reflectivity,  provided  data  is  available  over  the  portion  of  the  spectrum  where  activity  is 
significant . 

The  classical  models  discussed  have  limitations  which  arise  from  over-simplification  and  i 
improper  application  rather  than  lack  of  physical  credence.  The  multi-free/bound  electron 
model  as  shown  by  Roberts  has  been  successful  in  every  instance  where  a  comparison  with  data 
has  been  made.  This  model  suffers  from  its  cumbersome  form  and  large  number  of  parameters, 
giving  an  appearance  of  mathematical  expediency  rather  than  physical  reality.  Additional 
tests  of  this  model  require  optical  constants  data  or  complete  reflectivity  spectra.  The 
early  attempts  in  this  program  to  determine  the  model  parameters  from  an  incomplete  reflectiv- 
ity spectrum  can  now  be  judged  as  unfruitful  and  the  parameters  thus  obtained  are  mathematic- 
ally expedient  and  without  physical  significance. 

The  spectral  region  of  importance  to  thermal  radiation  applications  rarely  extends  to 
wavelengths  as  low  as  0.2p.,  since  only  0.18  percent  of  the  energy  from  a  blackbody  radiator  at 
3500  °K  is  at  shorter  wavelengths.  In  fact,  much  of  the  data  on  thermal  radiation  properties 
of  materials  do  not  extend  below  the  lower  limit  of  visible  light,  0.4  u.» 


40 


There  is  a  paradox  which  should  be  of  concern  to  workers  in  the  field.  For  practical 
applications,  data  are  required  from  the  visible  out  into  the  far  infrared,  while  for  theoret- 
ical computations,  data  are  required  over  a  wavelength  range  extending  into  the  extreme  vacuum 
ultraviolet.  Measurements  in  the  vacuum  ultraviolet  wavelength  range,  require  special 
experimental  techniques.  For  instance,  surface  contamination  may  exert  a  controlling  influence 
on  the  results  obtained.  It  is  not  possible  to  obtain  reliable  data  simply  by  extending  pres- 
ent techniques  used  at  the  longer  wavelengths. 

The  alternative  approach  to  specifying  optical  behavior  is  to  measure  two  parameters  at 
each  wavelength  of  interest,  as  was  previously  suggested.  This  would  require  different  or 
more  elaborate  techniques  than  are  now  in  use,  but  the  benefits  would  be  significant  and  last- 
ing. 

5.  Working  Standards  of  Normal  Spectral  Emittance 

5 . 1  Background 

Three  sets  of  working  standards  of  normal  spectral  emittance  have  been  prepared  and  cal- 
ibrated. Set  1  consisted  of  specimens  of  platinum-  13%   rhodium  alloy  and  oxidized  Inconel. 
JData  on  these  specimens  are  given  in  WADC  TR  59-510,  part  III,  C53l.  Set  2  consisted  of  a 
large  number  of  specimens  of  various  sizes  and  shapes  of  commerically  pure  platinum,  oxidized 
HLanthal,  and  oxidized  Inconel.  These  specimens  were  originally  prepared  for  distribution  by 
|ASD  to  Air  Force  contractors,  but  because  of  legal  questions,  they  were  returned  to  NBS  for 
sale  to  the  general  public*  Data  on  these  specimens  are  given  in  WADC  TR  59-510,  part  IV 
![54]«  Set  3,   like  Set  1,  consisted  of  specimens  of  platinum-  13$  rhodium  alloy  and  oxidized 
•llnconel.  They  were  prepared  for  use  by  RTD.  Data  on  these  specimens  are  given  in  NBS  Tech- 
nical Note  252  [55]. 

5.2  Comparison  of  Measured  Normal  Spectral 

Emittance  of  Platinum  to  Values  Predicted 
by  the  Hagen-Rubens  Equation 

The  data  referred  to  above  for  platinum  and  platinum-  13$  rhodium  alloy  working  standards 
were  compared  with  values  computed  from  the  resistivity  by  use  of  the  Hagen-Rubens  [19]  equa- 
tion (Drude-Zener) .  This  equation  relates  normal  spectral  emissivity,  e(0,\),  to  dc  resist- 
ivity, pnu  and  wavelength 

e(0A)  =  0.365(PmA)*  -  0.0464(pm/\)  (31) 

Equation  (31)  is  applicable  at  long  wavelengths  and/or  high  temperatures.  It  is  derived 
from  classical  electromagnetic  theory,  as  discussed  in  Section  4.4a  of  this  report. 

Figures  49  through  51  show  the  type  of  agreement  observed  for  platinum.  The  solid  lines 
Represent  the  average  measured  omittances  of  the  six  platinum  working  standards.  The  vertical 
lines  indicate  the  statistical  uncertainty,  expressed  as  the  total  standard  deviation,  due  to 
prror  of  measurement  and  difference  in  specimens.  The  two  broken  lines  represent  the  values 
bomputed  from  equation  (31),  using  resistivity  data  given  by  Caldwell  T36]. 

There  are  some  obvious  areas  of  mismatch  between  the  measured  and  computed  values.  The 
teximnm  near  9  microns  at  all  three  temperatures  could  be  caused  by  (a)  bound  electron  phenom- 

r,  (b)  an  instrumentation  error,  or  (c)  an  impurity  effect.  Recent  data  have  shown  that 
is  due  to  a  small  amount  of  silicon  impurity  in  the  platinum,  as  is  discussed  later. 

The  rise  in  emittance  values  at  wavelengths  near  15  microns  at  the  higher  temperature  is 
pelieved  to  be  an  instrumentation  error.  It  has  not  been  found  in  repeat  measurements  made 
bn  the  same  specimens  after  completing  certain  modification  in  the  spectrometer.  This  is 
pliscussed  in  Section  6.4. 

The  maxima  between  1  and  2  microns  have  previously  been  attributed  to  resonance  of  bound 
electrons,  although  data  reported  by  Roberts  [20]indicate  that  no  resonance  effect  is  to  be 
axpected  in  this  wavelength  region.  The  main  ma  may  be  due  to  instrumentation  error,  since 
they  occur  in  wavelength  regions  where  the  energy  is  low. 
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Figure  49.  Normal  Spectral  Emittance  of 
Platinum  Working  Standards  at  800  °K, 
Compared  to  Values  Computed  from  the 
Hagen-Rubens  Equation.  The  resistivity 
of  platinum  in  micro-ohm  centimeters 
is  29.1  at  800  °K. 
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Figure  50.  Normal  Spectral  Emittance  of 
Platinum  Working  Standards  at  1100  °K, 
Compared  to  Values  Computed  from  the 
Hagen-Rubens  Equation.  The  resistivity 
of  platinum  in  micro-ohm  centimeters 
is  38.1  at  1100  °K. 
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Figure  51.  Normal  Spectral  Emittance  of 
Platinum  Working  Standards  at  1400  °K, 
Compared  to  Values  Computed  from  the 
Hagen-Rubens  Equation.  The  resistivity 
of  platinum  in  micro-ohn  centimeters 
is  46.5  at  1400  °K. 
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Figure  52.  Normal  Spectral  Emittance  of 
Platinum-13^  Rhodium  Alloy  at  800  °K. 
The  resistivity  of  platinum-13^  rhodium 
alloy  in  micro-ohm  centimeters  is 
33.9  at, 800  °K. 
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Figure  53.  Normal  Spectral  Emittance  of 
Platinum  13$  Rhodium  Alloy  at  1100  °K, 
Compared  to  Values  Computed  from  the 
Hagen-Rubens  Equation.  The  resistivity 
of  platinum-13%  rhodium  alloy  in  micro- 
ohn  centimeters  is  41.9  at  1100  °K. 
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Figure  54.  Normal  Spectral  Emittance  of 
Platinum-13^  Rhodium  Alloy  at  1300  °K, 
Compared  to  Values  Computed  from  the 
Hagen-Rubens  Equation.  The  resistivity 
of  platinum-13/b  rhodium  alloy  in  micro- 
ohm  centimeters  is  46.7  at  1300  °K. 


At  1100  °K  (Figure  50)  a  reasonable  fit  of  the  experimental  data  at  wavelengths  beyond 
about  3.5  microns  was  obtained  with  the  values  computed  from  the  resistivity  of  platinum  at 
1100  °K  given  by  Caldwell  [36].  At  1400  °K,  however  (Figure  51),  a  resistivity  higher  than 
that  given  by  Caldwell  was  required  to  produce  even  an  approximate  fit,  and  at  800  °K,  a 
resistivity  lower  than  that  given  by  Caldwell  was  required.  These  results  indicate  that  equa- 

!  tion  (31)  cannot  be  used  with  published  resistivities  to  extrapolate  the  measured  omittance 

;  data  to  temperatures  above  1400  °K  or  below  800  °K. 

The  properties  of  the  platinum-  13$  rhodium  alloy  are  similar  to  those  of  pure  platinum. 
j  Normally,  the  presence  of  an  alloying  constituent  might  be  expected  to  introduce  changes  in 
the  properties.  However,  these  effects  are  minor  and  appear  as  a  slight  increase  in  omittance 
!  and  also  in  resistivity,  as  compared  to  pure  platinum.  The  experimental  data  are  compared  in 
I  Figures  52  through  54  with  values  computed  from  equation  (31) .  The  agreement  between  computed 
!  and  observed  values  is  about  as  good  as  was  the  case  for  platinum.  In  general,  however,  ap- 
j  plication  of  the  simplified  theory  to  alloys  should  be  done  with  caution,  since  an  alloy,  in 
I  most  cases,  can  no  longer  be  considered  as  a  simple  metal. 

The  observed  data  show  the  maximum  between  9  and  10  microns  that  was  observed  with  plat- 
jinum  at  each  temperature.  The  Tnaximum  between  1  and  2  microns  appears  only  in  the  platinum- 
1 1J%  rhodium  data  at  800  °K.  The  rise  in  omittance  at  wavelengths  near  15  microns  was  not 
I  observed  with  the  alloy  specimens. 

Several  attempts  were  made  by  spectrographic  analysis,  electron  diffraction,  and  X-ray 
diffraction  to  identify  any  impurity  that  might  be  present  in  the  platinum  and  platinum-  13$ 
rhodium  working  standards  of  normal  spectral  omittance  that  might  account  for  the  peak  in  the 
spectral  omittance  curve  near  9  microns.  Silica  was  suspected,  because  there  is  a  strong 
resonance  of  the  Si-0  bond  at  about  that  frequency.  However,  these  analyses  did  not  show 
the  presence  of  elements  other  than  platinum  or  rhodium  in  the  samples. 

Finally,  an  electron  probe  microanalysis  was  made  [ 37  ]  of  the  surface  of  a  platinum-13^ 
I  rhodium  working  standard,  which  indicated  the  presence  of  silicon  as  a  minor  constituent,  and 
1  a  trace  of  iron.  All  other  elements  of  atomic  number  above  12  were  found  to  be  absent. 

Area  X-ray  scanning  revealed  that  the  silicon  is  strongly  concentrated  at  the  grain 
boundaries,  as  shown  in  Figures  55>  56,  and  57,  and  is  absent  at  the  rest  of  the  surface. 
Strong  bluish  cathodoluminescence  observed  when  the  electron  beam  strikes  the  surface 
strongly  suggests  that  the  impurity  is  a  non-metallic  phase.  This  observation,  combined  with 
the  wavelength  of  the  emittance  peak,  indicates  that  the  impurity  is  probably  present  as  Si02» 
The  amount  of  the  impurity  present,  cannot  be  estimated  accurately  from  the  measurements,  but 
it  is  believed  to  be  certainly  less  than  0.1  percent,,  and  probably  less  than  0.01  percent. 

5.3  Comparison  of  Standards  Data  with  Literature  Values 

Since  the  optical  properties  of  a  metal  are  markedly  affected  by  the  conditions  of  the 
surface,  differences  in  values  reported  by  different  investigators  are  to  be  expected. 
Recently,  the  problem  of  surface  characterization  of  metals  [  38  ]  has  received  increasing 
attention.  The  complexity  of  the  problem  makes  it  extremely  difficult  to  characterize  com- 
pletely a  metal  specimen  so  that  data  obtained  in  one  laboratory  can  be  reproduced  on  other 
I specimens  of  the  same  material  prepared  and  measured  in  a  second  laboratory.  In  spite  of 
'these  difficulties,  comparison  of  data  obtained  on  the  working  standards  of  normal  spectral 
emittance  with  data  obtained  in  other  laboratories  on  similar  specimens  may  serve  a  useful 
purpose. 

Figure  58  shows  values  for  normal  spectral  emittance  of  platinum  at  1400  °K  obtained  in 
three  other  laboratories,  together  with  the  data  reported  for  the  working  standards.  The 
agreement  is  reasonably  good,  except  that  the  increase  in  emittance  at  14  to  15  microns  was 
not  reported  by  other  investigators  and  did  not  occur  in  more  recent  measurements  discussed 
in  Section  6.4. 

Figures  59,  60,  and  61  show  values  for  oxidized  Inconel  obtained  in  other  laboratories, 
together  with  the  data  on  the  working  standards.  In  general,  when  oxidation  times  were 
equivalent,  agreement  is  satisfactory. 
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Figure  55  •    Photomicrograph  of 
Platlnum-13^  Rhodium  Working 
Standard,  Unetched,  x  250. 


Figure  56.  Electron  Back- 
scattering  Image  of  Plat- 
inum-13^  Rhodium  Working 
Standard,  x  440. 


Figure  57.  Silicon  &*  X-ray 
Image  of  Platinum-13^ 
Rhodium  Working  Standard, 
x  440. 
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Figure  58.  Normal  Spectral  Quittance 
of  Platinum  Working  Standards  at 
1400  °K  Compared  to  Literature 
Values. 
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Figure  60.  Normal  Spectral  Emittance 
of  Oxidized  Inconel  Working  Standards 
at  1100  °K,  Compared  to  Literature 
Values. 
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Figure  59.  Normal  Spectral  Emittance 
of  Oxidized  Inconel  Working  Standards 
at  800  °K  Compared  to  Literature 
Values . 
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Figure  61.  Normal  Spectral  Emittance 
of  Oxidized  Inconel  Working  Standards 
at  1300  °K,  Compared  to  Literature 
Values . 


44 


1 

1                                     1                                    1 
4 

1      1      1      1 

— 

4 
4 

4 

4 

4 

4 

4 

4 

4 

— 

4 

<J              — 

4 

4 

4 

4 

4 

4 



4 

4          — 

4 

0 

4 

4 

4 

4 

4 

4                 

4 

0 

4 

<J 

4 

<J 

4 

<                             _ 

4 

4 

4 

4 

4 

4 



4 

4                             

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 



4 

4                 

4 

4 

4 

4 

4 

4                   

— 

4 

4 

4 

4 

4 

4 

4 

< 

4 

4 

— 

4 

4               — 

4 

4 

4 

4 

4 

4 

4 

4 

— 

— 

4 

4 

4fl 
4H 

v'  ~ 

- 

«g 

■  v- 

48 

*    > 

— 

<B 

0 

0 

0 

3 

I 

°°  0                _ 
0  0 

0 

_ 

0 
„0 

fl          - 

0 

_  0 

0 

0 

0 

0 

0 

0 

0 

c 

0 

0                                                 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

< 

1 

••> 

J>                                   — 

i 

8 

/ 

O 

°o 

0                                            — 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0              _ 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0                           — 

0 

0 

0 

0 

0 

0 

0 

0 

0 

„0 

0 

1 

1              1              1 

1       1       1       1 

o                      o 

10                   0                   0                  0 

tf>        ^           ro          esj           ~         0 

_  * 


P.-H 
CO   O 


001   X    30NV1JLIW3    nvai03dS   1VHMON 


001   X    N0I1VIA3Q    aavaNVis 


45 


A  brief  description  of  the  measurements  made  in  other  laboratories  is  given  in  Table  12. 

5.4  New  Working  Standards  of  Normal  Spectral  Emittance 

The  supply  of  working  standards  of  normal  spectral  emittance  for  sale  to  the  general 
public  became  seriously  depleted  by  mid-1964.  New  specimens  of  platinum-13^  rhodium  alloy, 
Kanthal,  and  Inconel  were  therefore  procured  and  prepared  for  calibration.  The  cost  of  pro- 
curement, preparation,  and  claibration  of  the  new  standards  is  being  financed  by  NBS  funds, 
but  because  of  the  great  interest  of  the  Air  Force  in  this  work,  and  its  close  connection  to 
other  work  on  this  project,  it  is  being  reported  here.  The  specimens  were  procured  in  the 
following  shapes  and  sizes: 

No.  of  Specimens  in: 


Size 

Inconel 

Kanthal 

Pt-13*  Rh. 

1/2"  discs 

20 

20 

20 

7/8"  discs 

22 

22 

22 

1"   discs 

20 

20 

20 

1  1/8"  discs 

16 

16 

16 

1  1/4"  discs 

16 

16 

16 

2"  x  2"  squares 

15 

15 

15 

1/4"  x  8"  strips 

7 

7 

7 

3/4"  x  10"  strips 

3 

3 

3 

1"  x  10"  strips 

3 

3 

3 

1"  x  1"  cylinders 

- 

- 

5 

The  platinum-13^  rhodium  specimens  were  0.035  inch  thick,  the  Kanthal  specimens  0.043 
inch  thick  and  the  Inconel  specimens  0.053  inch  thick.  The  platinum-13$  rhodium  cylinders 
are  to  be  calibrated  at  1400°,  1600°,  and  1800°  K,  by  the  rotating  cylinder  procedure.  The 
1/4  -  by  8-inch  strips  of  platinum-13^  rhodium  alloy  are  to  be  calibrated  at  800°,  1100°,  and 
1400 °K.  The  1/4-  by  8-inch  strips  of  Kanthal  and  Inconel  are  to  be  calibrated  at  800°,  1100 °,: 
and  1300°  K. 

The  platinum-13^  rhodium  specimens  were  washed  in  hot  tap  water,  then  in  distilled  water,' 
and  finally  in  ethyl  alcohol.  They  were  then  dried  in  air  and  placed  in  closed  containers, 
supported  by  the  ends  or  edges  only,  for  storage  prior  to  annealing. 

The  platinum-13^  rhodium  specimens  were  annealed  in  a  gas-fired  furnace.  The  strip 
specimens  were  hung  by  means  of  platinum-20$  rhodium  hooks  from  alumina  rods  in  the  furnace. 
The  disc  specimens  were  supported  by  the  edges  only  on  ceramic  forms  resting  on  a  rigid  cerami " 
slab.  The  square  specimens  were  stood  on  edge  on  a  flat  ceramic  form,  leaning  against  a  ce- 
ramic support.  All  of  the  specimens  were  finally  enclosed  in  a  ceramic  muffle.  The  tempera- 
ture of  the  furnace  was  raised  to  2830  °F  (1828  °K)  over  a  period  of  six  hours  and  held  at  tha 
temperature  for  one  hour,  after  which  the  fuel  was  cut  off  and  the  dampers  were  closed.  The 
specimens  were  allowed  to  cool  in  the  furnace  for  three  days.  They  were  then  removed  from 
the  furnace  and  stored  in  closed  containers,  supported  by  ends  or  edges  only. 

The  Inconel  and  Kanthal  specimens  were  roughened  by  sandblasting  with  60-mesh  fused 
alumina  grit  at  about  70  psi  air  pressure,  care  being  taken  to  obtain  a  uniformly  rough  sur- 
face on  all  specimens.  They  were  then  cleaned  ultrasonically  in  acetone,  rinsed  with  fresh 
acetone,  and  passivated  for  one  minute  in  10  percent  nitric  acid  at  109  °F  (316  °K),  rinsed 
in  distilled  water  and  then  in  freshly  distilled  acetone. 

The  Inconel  strip  specimens  were  suspended  from  a  metal  rack  by  means  of  platinum-20^ 
rhodium  hooks,  and  the  rack  was  inserted  into  a  cold  furnace.  The  disc  specimens  were  suppose 
by  the  edges  only  on  special  ceramic  supports  placed  on  the  hearth  of  the  furnace.  The  squart 
specimens  were  stood  on  edge  on  the  hearth  of  the  furnace,  leaning  against  a  ceramic  support.  . 
The  furnace,  heated  by  silicon  carbide  resistance  elements,  was  brought  to  a  temperature  of 
1952  °F  (1340  °K)  in  a  period  of  about  5  hours  and  held  at  that  temperature  for  24  hours, 
then  lowered  to  1520  °F  (1100  °K)  and  held  for  an  additional  24  hours,  after  which  the  power 
was  cut  off  and  the  specimens  were  allowed  to  cool  in  the  furnace  for  two  days.  They  were 
then  removed  from  the  furnace  and  stored  in  closed  containers,  supported  by  the  ends  or  edges 
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only.  The  Kanthal  specimens  were  oxidized  following  the  procedure  outlined  above  for  Inconel, 
except  that  the  heat  treatment  applied  for  was  400  hours  at  1952  °F  (1340  °K) .  The  specimens 
were  stored  in  a  similar  manner. 

Specimens  of  all  alloys  were  subjected  to  a  minimum  of  handling  after  the  initial  clean- 
ling.  When  handling  was  unavoidable,  metal  tongs  were  used,  or  rubber  surgical  gloves  were 
Iworn,  and  the  specimens  were  touched  only  on  the  ends  or  edges. 

6.  Normal  Spectral  Emittance 

6.1  Background 

Equipment  and  procedures  were  developed  for  the  measurement  of  normal  spectral  emittance 
,'of  specimens  of  metal  and  other  materials  that  could  be  heated  by  passing  a  current  through 
;them,  at  temperatures  in  the  range  of  800°  to  1400  °K,  and  at  wavelengths  in  the  range  of  1 
jto  15  microns,  as  described  in  WADC  TR  59-510,  part  IV,  1*54] .  During  the  preceding  contract 
year,  an  effort  was  made  to  extend  the  wavelength  range  out  to  35  microns  or  beyond.  Con- 
siderable progress  was  made  in  modifying  the  equipment,  but  satisfactory  operation  was  not 
achieved.  The  modifications  included  (l)  construction  of  an  air  purification  system  to  elim- 
inate or  reduce  atmospheric  absorption,  and  (2)  use  of  selective  filters  to  reduce  errors 
due  to  scattered  radiation. 

6.2  Modification  of  Equipment 

A  new  Reeder  vacuum  thermocouple  detector  was  procured  and  installed.  The  sensitivity 
of  the  new  detector  was  about  2.4  times  that  of  the  previous  detector.  The  entire  optical 
path  from  the  entrance  to  the  chopper  cover  to  the  detector  was  enclosed  in  a  light-tight 
and  air-tight  enclosure.  Provision  was  made  to  insert  or  remove  filters  from  the  optical  path 
without  opening  the  path  to  the  room  atmosphere.  It  was  found  there  was  significant  leakage 
of  room  air  into  the  recirculating  atmosphere  control  system  through  the  open  crankcase  of 
the  recirculating  compressor.  The  moisture  and  carbon  dioxide  entering  the  system  with  the 
room  air  soon  depleted  the  drying  towers  in  the  air  purification  system.  The  situation  was 
corrected  by  adding  a  laboratory  air  drier  between  the  compressor  and  the  air  purification 
train.  With  this  drier  in  operation,  the  useful  life  of  the  absorbent  in  the  drying  towers 
of  the  air  purification  system  has  been  increased  from  about  8  hours  to  over  two  weeks. 

6.3  Test  Procedure 

Preliminary  measurements  showed  that  the  techniques  developed  for  use  in  the  1-15  micron 
range  were  not  adequate  for  the  15-40  micron  range.  The  difficulties  in  obtaining  reliable 
data  in  this  long-wavelength  range  are  associated  with  (a)  low  energy,  (b)  scattered  radiation 
within  the  monochromator,  and  (c)  high  atmospheric  absorption.  All  of  these  difficulties 
need  to  be  corrected  or  compensated  for,  at  least  in  part,  if  reliable  data  are  to  be  obtained 
in  this  long-wavelength  range. 

At  1400  °K,  only  about  1  percent  of  the  total  energy  emitted  by  a  blackbody  radiator  is 
at  wavelengths  longer  than  16  microns.  At  lower  temperatures  the  amount  of  energy  in  this 
wavelength  range  is  less,  although  the  percentage  of  the  total  energy  is  greater.  The  reduc- 
tion in  energy  available  for  measurement  has  been  compensated  for  in  part  by  use  of  the  more 
sensitive  Reeder  vacuum  thermocouple  detector.  The  lower  available  signal  strength  resulting 
from  the  smaller  amount  of  energy  in  the  long-wavelength  range  requires  that  the  slit  servo- 
mechanism  be  adjusted  to  produce  a  constant  Io  signal  of  5  to  2.5  mv,  instead  of  the  10  rav 
normally  used  in  the  1-15  micron  range.  This  smaller  I  signal  results  in  some  reduction  in 
precision  of  the  measurements. 

The  scattered  flux  within  the  monochromator  is  largely  undispersed  flux  from  the  beam 
entering  the  monochromator,  and  the  major  portion  is  at  wavelengths  near  the  peak  of  the 
Planck  distribution  function  for  the  temperature  of  the  source.  It  can  be  eliminated  by  plac- 
ing a  selective  filter  over  the  entrance  aperture  to  the  monochromator.  Such  a  filter 
should  absorb  all  of  the  short-wavelength  radiation  (below  15  microns)  and  pass  all  of  the 
I  long-wavelength  radiation.'  An  alternative  approach  is  to  replace  the  diagonal  plane  mirror 
'near  the  exit  slit  of  the  monochromator  with  a  scatter  plate  that  reflects  the  long-wavelength 
I  radiation  specularly,  thus  focusing  it  on  the  detector,  while  reflecting  the  short-wavelength 
flux  diffusely,  so  that  most  of  it  is  lost  before  reaching  the  detector. 
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The  atmospheric  absorption  was  largely  eliminated  by  removing  from  the  optical  path  most 
of  the  carbon  dioxide  and  water  vapor. 

Two  techniques  were  developed  in  an  attempt  to  obtain  useful  data  in  the  15-40  micron 
wavelength  range.  In  the  first  technique  (A),  a  black  polyethylene  filter  was  placed  immedi- 
ately in  front  of  the  entrance  slit  to  the  monochromator.  This  filter  absorbs  essentially 
of  the  radiation  at  wavelengths  shorter  than  about  12  microns,  and  transmits  25  to  50  percent 
of  the  radiant  flux  in  the  15-40  micron  range.  In  addition,  the  slit  servomechanism  was 
adjusted  to  produce  an  I0  signal  of  2.5  mv  in  the  constant  I0  mode,  at  which  setting  a  constant 
I0  signal  was  produced  while  scanning  from  15  to  about  38.5  microns. 

The  second  technique  (B)  was  to  scan  the  15-38.5  micron  wavelength  range  in  three  por- 
tions: approximately  13-27,  27-32,  and  32-38.5  microns,  respectively.  In  the  first  portion, 
the  black  polyethylene  filter  and  a  normal  Io  setting  of  10  mv  were  used.  In  the  second  por- 
tion, the  filter  was  removed,  a  CaF2  reststrahlen  plate  was  inserted  in  place  of  the  plane  di- 
agonal  mirror,  and  the  Io  signal  was  adjusted  to  5  mv-  The  CaF£  reststrahlen  plate  is  a  good 
specular  reflector  for  radiation  at  wavelengths  longer  than  about  25  microns,  but  scatters 
radiation  of  shorter  wavelengths.  In  the  third  and  final  portion,  the  reststrahlen  plate  was 
used,  and  the  I0  signal  was  adjusted  to  2.5  mv.  It  was  thought  that  the  second  technique, 
while  somewhat  less  convenient  than  the  first,  should  give  greater  potentiometer  sensitivity 
in  the  first  two  portions,  and,  hence,  probably  better  precision,  but  it  did  not  approach 
expectations. 

The  test  procedures  were  evaluated  by  using  the  combination  of  a  50  percent  transmittance' 
sector-disc  attenuator  and  a  laboratory  blackbody  furnace  to  simulate  a  specimen  having  a  nor- 
mal spectral  emittance  of  0.50  at  all  wavelengths.  Four  tests  were  made  at  each  of  the  three 
temperatures,  800,  1100,  and  1300  °K,  using  each  of  the  two  procedures.  The  data  were  reduced 
at  23  wavelengths  in  the  12.4  to  26.8  micron  range,  15  wavelengths  in  the  26.4  to  32.6  micron 
range,  and  15  wavelengths  in  the  32  to  37.65  micron  range. 

The  data  in  Table  13  indicate  that  in  the  two  shorter  wavelength  intervals,  where  the  IQ 
signal  is  different  for  the  two  procedures,  the  over-all  scatter,  measured  as  the  average 
standard  deviation  at  the  23  and  15  wavelengths,  respectively,  was  smaller  for  procedure  A 
than  for  procedure  B  in  three  of  the  six  cases,  and  the  data  in  table  14  show  that  the  accuracy 
was  better  for  procedure  A  than  for  procedure  B  in  all  six  cases.  On  the  basis  of  these  re- 
sults, procedure  A  was  selected  for  use  in  future  measurements  in  the  15-38  micron  range. 

6.4  Measurements  on  Platinum 

The  normal  spectral  emittance  of  a  platinum  working  standard  was  measured  over  the  wave- 
length ranges  of  1-15  microns  and  13-28.5  microns,  at  1300  °K.  Three  separate  runs  were  made 
and  the  data  were  reduced  as  before.  The  average  emittance  values  and  standard  deviations  are 
plotted  in  Figure  62.  The  curve  for  the  average  emittance  values  shows  that  the  results  with 
the  sodium  chloride  prism  (1-15  M>)  and  those  with  the  cesium  bromide  prism  (12.8-28  n)  check 
well  within  the  known  precision  of  measurement  in  the  range  where  the  two  curves  overlap.  The 
data  on  the  standard  deviations  of  the  measured  values  indicate  that  the  precision  of  measure- 
ment is  about  the  same  in  the  two  ranges.  The  rise  in  emittance  in  the  13-15  micron  region 
that  was  found  in  the  original  calibration  (  Figure  58)  was  not  confirmed  in  the  measurements 
with  either  prism  (Figure  62). 

7.  Summary  and  Plans  For  Future  Work 

Construction  of  a  laser-source  integrating  sphere  reflectometer  was  completed.  The 
sphere  coated  with  3M  Velvet  white  paint  was  calibrated  for  integrating  and  reflecting  proper- 
ties at  0.632  microns,  and  preliminary  tests  indicated  that  the  shallow  cavity  technique  offer? 
advantages  for  the  calibration  of  integrating  sphere  ref lectometers .  A  heater  for  use  at 
temperatures  up  to  800  °K  was  constructed,  and  preliminary  tests  showed  that  emitted  flux  from 
the  hot  specimen  willnot  invalidate  reflectance  measurements  on  specimens  at  temperatures  up 
to  at  least  2500  °K,  pnd  perhaps  3000  °K.  Future  plans  include  (l)  modification  of  the  detect- 
or viewing  configu'-  -.on  to  reduce  errors  found  with  the  present  configuration,  (2)  recoating 
of  the  sphere  with  Darium  sulphate,  to  replace  the  3M  Velvet  white  paint,  (3)  calibration  of 
the  response  of  the  modified  reflectometer  at  0.632,  1.15  and  3.39  microns,  (4)  construction 
of  a  specimen  heater  to  operate  at  temperatures  up  to  1800  °K,  and  (5)  design  of  a  specimen 
heater  to  operate  at  temperature  up  to  3000  °K. 
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Calibration  of  the  reflecting  properties  of  the  ellipsoidal  mirror  reflectometer  is 
essentially  complete.  A  preliminary  error  analysis  indicates  that  the  instrument  is  capable 
of  measuring  absolute  reflectance  at  room  temperature  under  conditions  approximating  normal 
illumination  and  hemispherical  viewing  with  an  accuracy  of  at  least  2  percent,  and  perhaps 
1  percent  under  ideal  conditions.  Future  work  with  this  instrument  will  be  carried  on  under 
a  project  financed  by  the  National  Aeronautics  and  Space  Administration. 

A  brief  review  of  the  literature  on  the  relation  of  omittance  and  reflectance  to  other 
properties  of  metals  and  dielectrics  is  presented,  together  with  a  brief  summary  of  work  done 
in  predicting  the  properties  of  platinum  and  rhodium.  It  was  shown  that  the  complete  dis- 
cription  of  the  optical  behavior  of  any  material  requires  evaluation  of  two  parameters  at 
any  one  wavelength,  since  two  parameters  completely  characterize  the  form  of  the  electromag- 
| netic  wave  interacting  with  the  material.  As  an  alternative,  the  Kramer-Kronig  relation  pro- 
vides a  means  of  determining  optical  behavior  from  a  knowledge  of  one  parameter,  reflectivity, 
at  all  wavelengths  of  the  spectrum  where  optical  activity  is  significant.  No  further  work  is 
planned  on  the  relation  of  omittance  to  other  properties  of  materials. 

New  working  standards  of  normal  spectral  omittance  have  been  procured  on  NBS  funds  and 
prepared  for  calibration.  These  consist  of  platinum-13$  rhodium  alloy,  oxidized  Kanthal  and 
oxidized  Inconel.  The  work  of  calibrating  these  specimens  will  be  completed  in  the  early  part 
of  the  next  contract  year. 

The  normal  spectral  emittance  equipment  was  modified  to  operate  in  the  15  to  37.5  micron 
region.  This  completes  the  development  work  with  this  equipment. 
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9.  Appendix  A 
Kramers-Kronig  relation  for  phase  in  terms  of  reflectivity. 

If  V  -   _  '1  +  -,  is,  as  a  function  of  the  wavelength  X,  analytic  throughout  the  lower 

half -plane,  and  is  never  zero  there,  then  we  can  define  f (X)  =  IriV  to  be  analytic  in  the  lowei 
half -plane.  If  X0  is  any  positive  real  number,  let  C^  be  the  contour  illustrated,  where  the 
large  semicircle  has  radius  A  and  the  small  one  radius  l/A. 


Call  these  semicircles  C^  and  C^  respectively,  and  call  the  straight  part  of  the  contour 
C  ^.  By  the  Cauchy  integral  theorem 


r  h\)    .    r      r 

Jc.  *-\.     AG'  +  lc!'  + 


.!■• 


(A-l) 


If  f(X)  approaches  zero  uniformly  as  X  approaches  infinity,  it  is  clear  that  the  first  inte- 
gral goes  to  ^ero  as  A  increases.  The  second  integral  is,  for  large  A,  approximately 

f(kQ)  '  qii  \7\~  >      and  so  reduces  to  j^fCXo).  The  third  integral  approachee  the  negative  of 

the  Cauchy  principal  value  of  the  integral  of  f  W  from-°°to  °°,  which  we  denote  pf^  ±IM  d^ 


We  thus  have 


f<\>  =  +  £  rf^ 


(A-2) 


Now  if  we  write  V  =  veJ  ,  then  f(^)  =  IriV  =  lnv  +  j  6,  and  we  can  equate  imaginary  parts  in 
(A-2)  and  get 


\-\ 


(A-3) 


where  V^  is  the  reflectivity. 

By  a  simple  transformation,  we  can  obtain 


9(\  )  =  -  pT  ^lnv^)  -  lnv(-\)1  +  \o[lnv(\)  +  lnv(-X.)1 
o'   tt  -J  72       ~2 


*»  -  K 


and,  making  the  formal  assumption  that  v(  x)  =  v(-x), 


2\ 


(A-4) 


•cv--^»JA¥-> 


(A-5)  i 


\  -\ 
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ind  an  integration  by  parts  gives  us 

9(\  )   - 


1    p     v'(\) 

TT    J  0    V<\) 


.  \+\ 


In 


\-\ 


dX 


(A-6) 


a  form  which  is  sometimes  more  convenient. 

To  show  that  our  proposed  equations  for  the  reflectivity  satisfy  the  above  relations 
regardless  of  the  values  of  the  coefficients,  it  is  sufficient  to  show  that  the  V(X)  defined 
by  the  equations  is  analytic  in  the  lower  half  of  the  X  plane,  is  never  zero  there,  and  that 
for  real  X,  |  V(— X) I  =  |V(X)  I.  The  first  is  easy  to  show;  since  the  coefficients  are  all 
'positive,  it  is  easily  seen  that  each  term's  denominator  can  be  zero  only  for  X's  with  a 
Ipositive  imaginary  part.  The  last  condition  also  holds:  changing  the  sign  of  \,   when  \   is 
I  real,  has  the  sole  effect  of  changing  the  sign  of  the  imaginary  part  of  n2  and  so  of  nj  which 

leaves  HijiryJ  unaffected. 

Table  1.  Data  Obtained  with  3M  Velvet  yellow  paper  in 
First  Reflection  Cavity 


Viewing   , 
Trial    Configuration^/ 


L/R  Ratio2/ 


0.5 


0.75 


1.0 


1.5 


2.0 


Linearity 


■2/ 


Theoretical,  center^ 
point  illumination^' 

0.800 

0.640 

0.500 

0.308 

0.200 

0.751 

0.500 

0.253 

101 

Black  sleeve,  light 
pipe,  .030"  beam 

.832+1 

.-648+2 

.509±2 

.3H±1 

.191±1 

.745 

.493 

.257 

104 

Black  sleeve  only 

.832+5 

.683±1 

.561+3 

.332+2 

.189+1 

.749 

.498 

.251 

105 

White  sleeve,  opal 
glass  view  cap, 
light  pipe 

.822+4 

.563+0 

.426+0 

.272+0 

.178+0 

.747 

.498 

.251 

106 

Same  as  105,  but 
0.075"  beam 

.828+4 

.561+1 

•443±2 

.278+0 

.182+0 

.745 

.499 

.249 

107 

Long  black  sleeve, 
0.075"  beam 

.839±6 

.696+1 

.570+2 

.366+4 

.184±1 

— 

— 

.250 

108 

i"  baffle  colli- 
mator, 0.075"  beam 

.844 

.707 

.592 

.371 

.201 

— 

— 

— 

109 

k"   baffle  colli- 
mator, 0.075"  beam 

.842+2 

.702+1 

.578+2 

.392+1 

.173±1 

.746 

.496 

.247 

110 

Same  as  109,  EMI  5685 
tube  0.075"  beam 

.841 

.704 

.573 

.368 

.192 

.746 

.496 

.248 

111 

White  sleeve,  light 
pipe,  0.075"  beam 

.843±3 

.641±1 

.501+1 

.3U±1 

.200+1 

.745 

.496 

.250 

1/  See  figure  7- 

2/  +  values  indicate  total  spread  in  two  measurements.  Only  one  measurement  was  made  if  no 
+  value  is  given. 

3_/  Sector  disc  attenuator  check. 

hj    From  relation  f  =  l/[l  +  (L/R)2] . 
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Viewing 


Table  2.  Data  Obtained  with  Shallow  Cavity  Lined 
with  3M  Velvet  Paper 


L/R  Ratio 


Linearity^/ 


Trial  Configuration^  Lining  0.5      0.75    1.0     1.5     2.0 

102  Black  sleeve,    Tellow   .94211-2/  .86512  .80311  .67511  .56411  0.745  0.498  — 
light  pipe 

103  Black  sleeve    Tellow   .94711    .89712  .849+0  .69511  .56212   .746   .498 
only 

112   White  sleeve,    Yellow   .95013    .85011  .79313  .678+0  .56711   .745   .496 
light  pipe, 
0.075"  beam 

301    Same  as  112,    Green    .832+0    .65712  .542+3  -36611  .24511   -753   .491  0.248 

f3_/ 


302    Same  as  112 


Green    .81111    .60812  .47912  .30411  .19011   —     — 


1/  See  Table  1. 

2/  1  values  indicate  total  spread  in  two  measurements. 

3_/  First  reflection  cavity 

ij  Sector  disc  attenuator  check. 


Table  3. 

Absolute  Reflectance  Computed  1 
Shallow  Cavity  Measurements 

rom 

A. 

Method  1  Analysis,  F 

=  f . 

L/R 

Trial 

0.5 

0.75 

1.0 

1.5 

2.0 

102 
103 
112 
302 

.88013 
.89313 
.90017 
•  36213 

.85314 
.89713 
.830+2 
.146115 

.86011 
.903+0 
.850+3 
.270115 

.850+1 
.  863jl 
.853+0 
.31515 

.83311 
.83212 
.836+1 
.26515 

B. 

Method  2  Analysis,  f 

=  F. 

L/R 

Trials 

0.5 

0.75 

1.0 

1.5 

2.0 

101  -  102 
103  -  104 
111-112 
301  -  302 

.851+3 
.86713 
.874110 
.27414 

.84913 
.87913 
.83012 

.333113 

.85612 
.884+0 
.84913 
.355113 

.84811 
.85611 
.850+0 
.329+5 

.83411 

.e: 

.83511 
.30517 

C. 

Method  3  Analysis,  F 

for  first  reflection,  f  for  all 
L/R 

others . 

Trials 

0.5 

0.75 

1.0 

1.5 

2.0 

101  -  102 
103  -  104 
111-112 
301  -  302 

.78612 
.80212 
.80519 
.260+0 

.842+3 

.837+3 

.82912 
.356110 

.85111 
.845+0 
.84912 
.365111 

.84711 
.84- 
.847+0 
.32913 

.840+1 

.84:z- 

•8"~z- 
.3: 
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Table  4.  Reflectance  of  Various  Materials  at  0.6328  Micron, 
Relative  to  3M  Velvet  Yellow  Paper 


Specimen 


Relative  Reflectance 


Deviatio 


ni/ 


3M  Velvet  yellow  over 
white  bond  (standard) 

3M  Velvet  yellow  over 

brass 
3M  Velvet  green  over 
white  bond 

3M  Velvet  white  over 
white  bond 

3M  Velvet  black  over 
white  bond 

ZrOfe   #1 

Zr02  #2 

CaZr03  #1 

CaZrOfl  #2 

Th02  #1 

Th08  #2 

Al^Oa   #1 

A1S03  #2 

ZrOa  #1 

Zr02  #2 

Rhodium  Mirror  » 

Aluminum  Mirror 

Gold  Mirror 

Silver  Mirror 

Porcelain  Enamel  #1 
#2 
#3 

#4 

#5 
#6 

#7 
#8 


Laser 
Reflectometer 

1.000 


0-373±l«^ 

1.02516 

0.044±1 
1.040+2 
1.039±2 
0.987±2 
0.953+3 
1.065+5 
1.033±2 
1.040+1 

1.039±2 
1.063±1 
1.067±1 
0.979±1 
1.119±1 
1.166+4 
1.208+5 
0.941+0 
0.081+1 
0.395±2 
0.092+1 
1.075±3 
0.187+1 

0.838+1 
0.287+1 


Commercial   , 
Ref le  ct  omet  er-^ 

1.000 

0.991 

0.364 

1.030 

0.035 

1.039±1 
1.060+1 

0.992+2 

0.972+4 

1.079±3 

1.050+2 

1.020+30 

1.028+30 

1.060+2 

1.060+3 

0.955±26 


0.915 
0.064 
0.373 

1.072 

0.167 

0.825 
0.263 


Total 


0.009 

2.4 

0.005 

-0.4 

0.009 

25.7 

0.001 

0.1 

0.021 

-2.0 

0.005 

-0.5 

0.019 

-1.9 

0.014 

-1.2 

0.017 

-1.6 

0.020 

1.9 

0.011 

1.1 

0.003 

0.2 

0.007 

0.6 

0.024 

0.026 

2.8 

0.017 

35.9 

0.022 

5.9 

0.003 

0.3 

0.020 

10.7 

0.013 

1.58 

0.024 

9.13 

1/  Relative  to  Commercial  Reflectometer  Value 

2/  Data  without  +  value  is  from  G.E.  Recording  Spectrophotometer.  The  +  value  is 
total  spread  of  two  readings  of  Cary  Model  14  Spectrophotometer. 

3_/  The  +  value  is  total  range  of  three  measurements. 

4/  Not  the  same  specimen. 
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Table  5.  Parameters  and  Relations  Descriptive  of  Wave-  Matter  Interaction  [21] 
e*  =  e'  -  Je" 


Complex  Permittivity 
Complex  Permeability 


i*  = 


MO 


Complex  Propagation  Factor  y   =  a   +  J  3 


Complex  Dielectric  Constant  K*  =  —  =  k'  -j  K1 


:0 


Complex  Index  of  Refraction  n*  =  n  -  njk 


Intrinsic  Impedance         z*  =  V  ^ 


Y  -  J  co  (e*  jx*)5 


n*  =  VK* 


e'  -  Real  Permittivity 
e"  -  Loss  Factor 

Valid  Assumption  for  Almost  All  Applica- 
tions 

a  -  Attenuation  Factor 
13  -  Phase  Factor 

i 
K  -  Relative  Permittivity  or  Real  Die- 
lectric Constant 

ti 
K  -  Relative  Loss  Factor 

n  -  Index  of  Refraction,  Ratio  of  Phase  j 

Velocities 
k  -  Index  of  Absorption 

Z  -  Ratio  of  Electric  to  Magnetic  Field 
Vectors 

E*  =  eo(l-jZ) 


Table  6.  Typical  Values  of  Parameters  Descriptive  of  Wave-Matter  Interaction  f27] 


Material, 
Temperature 
and  Wavelength 
Theoretical 
Conditions 

Complex 
Propagation 
Factor  y 

Complex 
Dielectric 
Constant  K 

Complex 
Index  of 
Refraction  n 

Intrinsic 
Imped ence  Z 

Normal 

Spectral 

Emissivity 

Q 

(n-1) 

pGt1) 

K' 

K" 

n 

k 

z'fa) 

Z"(fi) 

eQ,\ 

Copper 

300  K,0.5  n 

35 

13 

-6.69 

5.74 

1.03 

2.69 

389 

1048 

•  45 

Tungsten 
300  K,0.5  ll 

40 

44 

2 

22 

3.5 

0.91 

1300 

1190 

•  50 

Nickel 

300  K,0.5  m> 

40 

20 

-7-5 

10 

1.6 

2.0 

603 

1206 

•  40 

Mercury 

300  K,0.589  m, 

47 

29 

-16.8 

14.3 

1.62 

2.72 

610 

1660 

.25 

Synthetic 

Sapphire 

300  K,4  \l           2.10  x 

2.63 

2.81 

4-49  x 

1.675 

8.00  x 

630 

5.03  x 

10~6 

10-6 

10-7 

10-4 
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Table  7-  Classical  Models  for  the  Prediction 
of  the  Bnissivity  of  Metals 


e0,x  =  -365  (PmA)2  "  -°^  (PinA)  +  •  •  ■ 


Also  called  the  Drude- 
Zener  relations;  appli- 
cable at  long  wavelengths 
and/or  high  temperature; 
derived  from  classical 
electromagnetic  theory 
assuming  negligible 
phase  change  of  electric 
field  within  the  metal. 


e0,t  =  -5736  (pmT)5  -.1769  (pmT) 


Obtained  by  integration 
of  Hagen-Rubens'  first 
term  with  Planck  distri- 
bution curve,  when  com- 
bined with  p  vs  T  linear- 
ity relation  then  radiant 
energy  is  proportional  to 
the  fifth  power  of  the 
temperature. 


^0,t 


•  5736  (pmT)2  -.1769  (pmT)  +  .058  (p^)3/2   Obtained  by  integration 

of  Hagen-Rubens1  first 
and  second  terms  with 
Planck  Distribution  curve. 


:h,t  =  -751  (PmT)  +  '  • PmT  <•! 


Using  the  Fresnel  Relation 
for  angular  dependency, 
hemispherical  emissivity 
is  obtained. 


*h,X 


=    .476   (p   /\)2  -.148   (p   /X)   (Xv  /\<.5 


\* 


.442    (p   A)*  -.0995    (P   /\)   ,5<-P   A<2.5 
m  m        >  m 


h.t 


=   .751    (p  YV  -.396    (p  T)   o<P  T  <.2 

m  m     '  m 


With  graphical  integration 
techniques,  the  limit  of 
applicability  of  the 
Davisson-Weeks  Model  is 
entended. 


*,* 


=  .698  (p  T)"8  -.266  (p  T)t   .2  <  pT  <.5 

m  m  > 


where 


'0.\ 


25,  (-D n+i «rpmy2  PT«7.8 

n=l  rr.-^    >      m    - — — 


6  =  [(1  +  y2)Vl%,  y  = 


30\ 

2TTcm 

\p   Ne2' 
m 


i  +  r 


m 
m  ~  Ne2T 


Assumed  the  Drude  Theory 
with  finite  relaxation 
time;  expressions  for 

e0,t>  eh,X>  ***■  eh,t  are 
also  developed. 
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Table  8.  Theoretical  Models  for  Prediction 
of  the  Optical  Properties  of  Metals 


Free-Electron 


Assumes  free  electrons  of  one  type  are  subjected  to  sinusoidal  electric  field  and 
damped  by  a  viscous  force. 


(±- 


-\ 


K  =  1 


(h2\     1  +  lUJ 

V) 


V 


1  + 


<£? 


Multi-Free/Bound 

Assumes  several  types  of  free  and  bound  electrons  are  subjected  to  viscous  damping  arte 
elastic  restoring  forces. 

2 


N  \~     \~                            \2   q2 
K*  =  1  +     °      E  **     0h  -  — 2o 


2   2 

bh  Qh 


E  fg_fg_ 


4tt  m«0c 


^-4-£w 


2   8   \c   "  J*- 
4tt  m«0c   B    fg   J 


Roberts 


Same  form  as  above,  but  no  significance  can  be  attached  to  constants.  Shown  to  be 
valid  for  *n  metals  and  temperatures  considered. 


K  =  1  +  E  -= r 

m   ,  *■   -  _  *■ 


K  \ 

an 


\   +  j  6  X  \ 

sm      m  sm 


2nce, 


a 

E  "-^T 

n  X.   -  JA. 
rn 


Table  9-  Comparison  of  the  Roberts'  Multi- 
Free/Bound  Electron  Model  Constants  [27] 


Electron  Type 


Free 


Roberts ' 
Constants 

on 


Multi-Free/Bound 
Parameters 


A 


It 


2tt  c  m  A>, 


Free 


Bound 


Bound 


Bound 


rn 


om 


sm 


sm 


l* 


q% 


hnjQfc_ 


/  2     2 
4tt  me   c 


Lsm 


XOlAbh 


Physical  Significance 


Contribution  to  conductivity  by  free 
electrons,  type  g,  when  X-.»,  the  summa- 
tion of  contributions  by  all  types  gives 
DC  electrical  conductivity. 

Relaxation  wavelength  of  free  electorns 
giving  an  indication  of  time  between 
electron  collisions. 

Proportional  to  the  ratio  of  Coulomb 
repulsive  force  between  bound  electrons 
and  the  elastic  restoring  force. 


Resonance  wavelength  of  bound  electrons 
without  damping  forces  present. 

Ratio  of  resonance  wavelength  to  relax- 
ation wavelength  under  condition  where 
damping  and  restoring  forces  are  respec- 
tively absent. 
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Table  12  Brief  Description  of  Methods  of  Measurement  and 
Sample  Preparation  for  Cited  Literature  Values 


iterial  Authors 

3latinum    Schatz  &  McCandless  [41] 


Platinum     Maki,  Stair  &  Johnson  [40] 


Oxidized 
Inconel 


Oxidized 
Inconel 


Maki,  Stair  &  Johnson  [40] 


Blau,  et  al.  [42] 


Platinum     Price  [39] 


Oxidized 
Inconel 


Slemp  &  Wade  [44] 


Method/Apparatus 

Double-beam  recording 
spectrophotometer  with 
transfer  optics  to 
view  separate  specimen 
and  blackbody  furnaces. 

Single-beam  recording 
spectrophotometer  with 
transfer  optics  and 
traveling  optical 
bench  to  alternately 
view  blackbody  furnace 
and  resistance  heated 
specimen. 

Same  as  above. 


Single-beam  recording 
spectrophotometer  with 
specimen  imbedded  in 
wall  of  hollow  silicon 
carbide  heating  element 
with  cavity  to  serve 
as  blackbody. 

Single-beam  recording 
spectrophotometer  with 
cylindrical  specimens 
electrically  heated  in 
vacuum  with  lateral 
slit  serving  as 
blackbody. 

Single-beam  recording 
spectrophotometer  views 
specimen  as  it  rotates 
past  viewing  port  of 
blackbody  cavity 
furnace. 


Specimen  Description 

Commercial  foil, 
no  additional 
description  given. 


Commercial  grade 
strip  (85  x  ;  x  |") 
annealed  at  1525  °K 
for  one  hour. 


Commercial  Inconel 
strip,  sandblasted 
and  heated  at  ele- 
vated temperature 
for  six  hours. 

No  details  given 
other  than  oxidized 
in  air. 


Commercial  99.8$ 
purity  sheet. 


Oxidized  in  air 
for  20  minutes. 
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Temperature 


Table  13 .   Table  of  Average  Standard  Deviations 


Wavelength  Range ,   u 


12.40  - 

■  27.30 

26.40 

-  32.40 

32.0  - 

37.65 

Procedure 

Pro. 

:edure 

Procedure 

A 

B 

A 

B 

A 

B 

800  °K 

0.431 

0.603 

0.325 

0.250 

0.586 

0.564 

1100  °K 

.398 

•  345 

.326 

•  923 

•  593 

.441 

1300  °K 

•  402 

•  334 

.388 

.504 

•  467 

•  351 

Tempe'rature 


Table  U+.     Table  of  Average  Measurer '•nts 
Wavelength  Range,  p, 


12.40  - 

-  27.30 

26.40  - 

32.40 

32.0  - 

37-65 

Procedure 

Procedure 

Procedure 

A 

B 

A 

B 

A 

B 

800  °K 

52.00 

52.23 

51.68 

52.43 

52.21 

53-27 

1100  °K 

50.53 

51-78 

50.52 

53-21 

50.66 

53-44 

1300  °K 

50.56 

50.66 

50.58 

51.18 

50.92 

51-56 
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